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ABSTRACT 

Raman scattering by H2 in Neptune’s atmosphere has significant effects on its reflectivity for A < 0.5 
/im, producing baseline decreases of ~ 20% in a clear atmosphere and ~ 10% in a hazy atmosphere. 
However, few accurate Raman calculations are carried out because of their complexity and compu¬ 
tational costs. Here we present the first radiation transfer algorithm that includes both polarization 
and Raman scattering and facilitates computation of spatially resolved spectra. New calculations 
show that Cochran and Trafton’s (1978, Astrophys. J. 219 , 756-762) suggestion that light reflected in 
the deep CH4 bands is mainly Raman scattered is not valid for current estimates of the CH4vertical 
distribution, which implies only a 4% Raman contribution. Comparisons with lUE, HST, and ground- 
based observations confirm that high altitude haze absorption is reducing Neptune’s geometric albedo 
by ^6% in the 0.22-0.26 pm range and by ^13% in the 0.35-0.45 pm range. A sample haze model 
with 0.2 optical depths of 0.2-p,m radius particles between O.I and 0.8 bars fits reasonably well, but 
is not a unique solution. We used accurate calculations to evaluate several approximations of Raman 
scattering. The Karkoschka (1994, Icarus 111 , 174-192) method of applying Raman corrections to 
calculated spectra and removing Raman effects from observed spectra is shown to have limited appli¬ 
cability and to undercorrect the depths of weak CH4 absorption bands. The relatively large Q-branch 
contribution observed by Karkoschka is shown to be consistent with current estimates of Raman cross 
sections. The Wallace (1972, Astrophys. J. 176 , 249-257) approximation, produces geometric albedo 
values ~5% low as originally proposed, but can be made much more accurate by including a scattering 
contribution from the vibrational transition. The original Pollack et al. (1986, Icarus 65 , 442-466) ap¬ 
proximation is inaccurate and unstable, but can be greatly improved by several simple modifications. 

A new approximation based on spectral tuning of the effective molecular single scattering albedo 
provides low errors for zenith angles below 70°in a clear atmosphere, although intermediate clouds 
present problems at longer wavelengths. 

Subject headings: Neptune, Neptune Atmosphere, Spectrophotometry, Radiative Transfer 


1. INTRODUCTION 

Because Neptune’s atmosphere has a relatively low 
burden of aerosols, its reflected spectrum is strongly in¬ 
fluenced by both Rayleigh scattering and Raman scatter¬ 
ing by molecular hydrogen. Rayleigh scattering induces 
polarization that can significantly modify the reflected 
intensity (Mishchenko et al. 1994), accurate computa¬ 
tion of which presents the very large burden of solv¬ 
ing the vector radiation transfer equation. Sromovsky 
(2004) discusses that problem and a new approximation 
method applicable to low phase angles. Accurate treat¬ 
ment of Raman scattering is also a computational burden 
because photons incident at one wavelength lose some en¬ 
ergy to rotating and/or vibrating the hydrogen molecule 
and reappear at longer wavelengths. Computation of re¬ 
flectivity at one wavelength thus requires accounting for 
contributions from Raman scattering at shorter wave¬ 
lengths. In addition, because the Raman source function 
varies continuously with optical depth, otherwise homo¬ 
geneous layers become inhomogeneous, requiring many 
more layers to achieve an accurate characterization of 
the atmosphere. 

To avoid the computational burden of rigorous Ra¬ 
man scattering calculations, several different approxi¬ 
mations have been employed. Baines and Smith (1994) 
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used an approximation suggested by Wallace (1972), in 
which the Raman cross sections for rotational transi¬ 
tions are treated as conservative scattering because the 
wavelength shifts are relatively small, while the cross 
section for the vibrational transition, which involves a 
much larger wavelength shift, is treated as an absorp¬ 
tion. This approximation does not produce the sharp 
spectral features characteristic of Raman scattering and 
is of uncertain accuracy. Pollack et al. (1986) used an al¬ 
ternate approximation in which the Rayleigh scattering 
cross section at a given wavelength is scaled by the solar 
irradiance ratio at the shifted and unshifted wavelengths. 
This approximation does produce Raman spectral fea¬ 
tures, but the accuracy is not well known, and can cre¬ 
ate conservation problems by allowing single-scattering 
albedo values exceeding unity. Karkoschka (1994) pre¬ 
sented a method for correcting observations to remove 
Raman scattering and for converting calculations that 
ignored Raman scattering to spectra that approximately 
matched spectra that included Raman scattering. That 
method can add or remove Raman spectral features and 
was applied to observations of Saturn, Jupiter, Uranus, 
and Neptune, but was never tested for accuracy or gen¬ 
erality. 

A number of calculations of Neptune’s geometric 
albedo have been made that do account for the ba¬ 
sic physics of Raman scattering. Cochran and Trafton 
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(1978) implemented an iterative algorithm in which Ra¬ 
man scattering is first treated as an absorption. After 
solving the scalar radiation transfer equation at each fre¬ 
quency grid point, the photon loss is computed from the 
radiation field. On the next iteration the lost photons are 
added back as source terms at the shifted wavelengths 
appropriate to each Raman transition. They achieved 
convergence after three iterations. Their mainly low res¬ 
olution results are of limited utility however, because the 
atmospheric structure they assumed is so different from 
our current understanding. Their claim that the resid¬ 
ual intensity in the cores of the strong methane bands 
could be entirely explained by Raman scattering will be 
shown to be invalid because of their assumed CH 4 mixing 
ratio profile. The first model calculations displaying ex¬ 
tensively detailed Raman spectral features in Neptune’s 
atmosphere are those of Courtin (1999), who made use 
of the two-stream code of Toon et al. (1989) to speed 
the solution of the radiative transfer equation. But this 
method can deviate from exact solutions by 10-15%, does 
not account for polarization, and is not usable for study¬ 
ing center-to-limb variations. A more rigorous method 
was used by Betremieux and Yelle (1999), based on the 
DISORT radiative transfer code (Stamnes et al. 1988). 
But they presented only results for Jupiter and omitted 
polarization effects. 

This paper presents a new method for accurate com¬ 
putation of Raman scattering that includes polarization 
in the context of Neptune’s atmosphere and makes pre¬ 
liminary applications of that method to resolve several 
significant issues. The next section reviews the basic 
physics of Raman scattering. That is followed by a dis¬ 
cussion of methods for accurate computation of Raman 
scattering. Sample computations are then presented to 
characterize the basic features of Raman scattering on 
Neptune. Comparisons are made with various prior cal¬ 
culations, and with HST and groundbased observations, 
to assess the degree of haze absorption required in Nep¬ 
tune’s atmosphere. The final section evaluates past ap¬ 
proximations, discusses how they can be generalized and 
improved, and presents the new approximation and its 
performance. 

2. THE PHYSICS OF RAMAN SCATTERING 
2 . 1 . Hydrogen Energy Levels and Transitions. 

The hydrogen molecule can exist in two nuclear spin 
states. Ortho states have parallel nuclear spins and odd 
total angular momentum quantum numbers (J = 1, 3,...) 
with a degeneracy of 3(2 J -|- 1). The para states have 
antiparallel nuclear spins and even angular momentum 
quantum numbers (J = 0,2,..) with a degeneracy of 
(2J -I- 1). The equilibrium population of these states 
follows the Boltzmann distribution, so that the fraction 
of molecules with angular momentum J is given by 

00 

Peq{J) = d{J) exp{-E{J)/kT)/ '^{d{J) exp{-E{J)/kT)) 

j=o 

( 1 ) 

where d(J) is the degeneracy, E{J) is the energy above 
the ground state, k is the Boltzmann constant, and T 
is absolute temperature. Using published expressions 
for the energy levels (Farkas 1935; Massie and Hunten 
1982) we obtain the fractional populations given in Ta- 


TABLE 1 

Fractional Populations of H2 Rotational States for 
Equilibrium H2. 


T(K) 

0 

II 

J=1 

II 

to 

J=3 

II 

fp 

fo 

50 

0.7704 

0.2294 

0.0001 

0.0000 

0.0000 

0.771 

0.229 

75 

0.5173 

0.4798 

0.0029 

0.0000 

0.0000 

0.520 

0.480 

100 

0.3747 

0.6135 

0.0115 

0.0003 

0.0000 

0.386 

0.614 

125 

0.2947 

0.6784 

0.0250 

0.0018 

0.0000 

0.320 

0.680 

150 

0.2450 

0.7080 

0.0410 

0.0059 

0.0000 

0.286 

0.714 

175 

0.2112 

0.7178 

0.0574 

0.0135 

0.0001 

0.269 

0.731 

200 

0.1865 

0.7157 

0.0729 

0.0245 

0.0004 

0.260 

0.740 

225 

0.1673 

0.7061 

0.0869 

0.0387 

0.0009 

0.255 

0.745 

250 

0.1520 

0.6919 

0.0990 

0.0552 

0.0017 

0.253 

0.747 

275 

0.1394 

0.6749 

0.1092 

0.0732 

0.0028 

0.251 

0.749 

300 

0.1287 

0.6564 

0.1177 

0.0919 

0.0043 

0.251 

0.749 


Note: fp and denote fpara and fortho respectively. 


TABLE 2 

Fractional Populations of H2 Rotational States for 
Normal H2. 


T(K) 

0 

II 

•—5 

J=1 

C-l 

II 

to 

J=3 

J=4 

fp 

fo 

50 

0.2500 

0.7500 

0.0000 

0.0000 

0.0000 

0.25 

0.75 

75 

0.2486 

0.7500 

0.0014 

0.0000 

0.0000 

0.25 

0.75 

100 

0.2426 

0.7496 

0.0074 

0.0004 

0.0000 

0.25 

0.75 

125 

0.2305 

0.7480 

0.0195 

0.0020 

0.0000 

0.25 

0.75 

150 

0.2142 

0.7438 

0.0358 

0.0062 

0.0000 

0.25 

0.75 

175 

0.1965 

0.7362 

0.0534 

0.0138 

0.0001 

0.25 

0.75 

200 

0.1795 

0.7251 

0.0702 

0.0249 

0.0004 

0.25 

0.75 

225 

0.1640 

0.7109 

0.0851 

0.0390 

0.0008 

0.25 

0.75 

250 

0.1504 

0.6944 

0.0979 

0.0554 

0.0016 

0.25 

0.75 

275 

0.1386 

0.6762 

0.1086 

0.0733 

0.0028 

0.25 

0.75 

300 

0.1283 

0.6570 

0.1174 

0.0920 

0.0043 

0.25 

0.75 


Note: ip and fo denote ipara and fortho respectively. 


ble 1 , where fpara is the total fraction of molecules with 
even J and fortho is the total fraction of molecules with 
odd J. This assumes that ortho and para states can ex¬ 
change energy. However, the time scale for equilibration 
by means of bimolecular collisions is of the order of years 
(Massie and Hunten 1982) and radiative energy exchange 
is forbidden by selection rules. Thus it is also meaning¬ 
ful to consider a different kind of equilibrium condition 
in which ortho states equilibrate separately from para 
states. This is relevant when hydrogen equilibrates at 
high temperature, then is lifted by convection to higher 
altitude where the fraction of para and ortho molecules 
remain fixed (over short time scales) but equilibration 
of energy levels within each sub-population does take 
place. For “normal” H 2 , which is defined by the high 
temperature equilibrium value of fpara = 0.25, the sub¬ 
population distributions are given by 

p in-/ exp{-E{J)/kT)/Zp, J even 
norm[J) \3d{J)exp{-E{J)/kT)/Zo, J odd ^ ’ 

where the partition functions Zp and Zq are summations 
of d{J) exp{—E{J)/kT) carried out over even J and odd 
J respectively. This distribution is given in Table 2 as a 
function of temperature. 

The populations for both equilibrium and normal hy¬ 
drogen are plotted in Fig. [U both as a function of tem¬ 
perature and as a function of pressure in Neptune’s at¬ 
mosphere. At T > 300 K, fpara approaches 0.25 and 
the ortho/para ratio approaches the 3/1 ratio expected 
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from the nuclear spin degeneracy. But at the low tem¬ 
peratures in the upper troposphere and stratosphere of 
Neptune fpara can be much larger, and for pressures less 
than a few bars, only J = 0 and J = 1 ground states 
need to be considered. 




Population Fraction for Neptune 


Fig. 1. — TOP: Fractional populations of angular momentum 
states as a function of temperature for equilibrium H 2 (thick lines) 
and for normal H 2 (thin lines). BOTTOM: Population fractions 
vs pressure in Neptune’s atmosphere. 

Radiative transitions must satisfy selection rules AJ = 
—2,0, +2 (Hollas 1992); the corresponding transitions are 
named 0{J), Q{J), and S{J), where J is the angular mo¬ 
mentum quantum number of the initial state. The O and 
S branches can involve changes in rotational energy and 
vibrational energy, while the Q branch involves changes 
in only vibrational energy. For the O branch interac¬ 
tions the scattered photons have more energy than the 
incident photons. But for the O branch to play a signif¬ 
icant role, there needs to be a significant population in 
rotational states with J > 2, which is not the case for 
the upper troposphere of Neptune. Thus, we here ignore 
the O branch. 

2 .2. Raman Scattering Cross Sections 

The wavelength-dependent Raman and Rayleigh cross 
sections per molecule (Fig. [5]) are computed using fits 
given by Ford and Browne (1973). Cross sections at 0.4 
/rm are given in Table 3, along with transition energy 
energies expressed as corresponding wavenumber shifts. 
These cross sections for the S and Q Raman transi- 


TABLE 3 

Raman transitions, wavenumber shifts, and 0.4-/xm cross 

SECTIONS PER MOLECULE. 






OA-fim Cross 

Transition( J) 

AJ 

Av 

Aiy, cm ^ 

section, cm^ 

R(0) 

0 

0 

354.69 

3.575x10-^^ 

R(l) 

0 

0 

587.07 

3.635x10-2'^ 

S{0) 

-f2 

0 

354.69 

1.104x10-28 

S(l) 

+2 

0 

587.07 

0.642x10-28 

Q(0) 

0 

1 

4162.06 

0.344x10-28 

Q(1) 

0 

1 

4156.15 

0.369x10-28 

Si(0) 

-f2 

1 

4498.75 

0.070x10-28 

Si{l) 

-f2 

1 

4713.83 

0.041x10-28 

Q 



4161.00 

0.412x10-28 


NOTES: The R{J) transitions refer to Rayleigh scattering 
cross sections, which do differ slightly with ground state 
quantum number. The subscript in Si(J) refers to the 
vibrational quantum number change. The unsubscripted 
Q refers to the sum of the four prior cross-sections divided 
by two (the last term in Eq. [Sjl . 

tions are 46-68% larger than those given by Cochran and 
Trafton (1978). 



Fig. 2.— TOP: Raman cross sections/molecule vs wavelength 
for J = 0 ground states (thick lines) and for J = 1 ground states 
(thin lines), where refers to the AJ = 2, Au = 1 vibrational 
transition. BOTTOM: Ratio of Raman to Rayleigh cross sections. 

Note in Table 3 that the <5(0) -I- 5'i(0) cross section 
sum is only 1% larger than Q[\) + S'i(l) sum at 0.4 ^m. 
Fig. [5] shows that this close match holds true at all wave¬ 
lengths of interest. These transitions also all have similar 
wavenumber shifts. In fact, the shifts of Q(0) and Q(l), 
which have the largest cross sections, differ by only 6 
wavenumbers, which is several times smaller than the 
spectral resolution that we need to model existing obser¬ 
vations. For these transitions we can thus follow Courtin 
(1999) by ignoring fpara and simply taking the effective 
cross section as the average of J = 0 and J = 1 transi- 
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tions and using a single wavenumber shift of 4161 cm“^. 
The inclusion of the terms correctly accounts for the 
absorption of photons at the incident wavelength, but 
does not transfer the scattered photons to exactly the 
correct longer wavelength. However, this results in a 
negligible error. 

The effective scattering cross section per molecule for 
each transition depends on the fractional population of 
the initial state as well as the transition cross section. For 
Neptune’s the upper troposphere we can make the ap¬ 
proximation that there are only two ground states: J = 0 
and J = 1. This means that fpara can be taken as the 
fraction of J = 0 molecules, and fortho as the fraction 
of J = 1 molecules. We can then write the total Raman 
scattering cross section as 


3. RADIATION TRANSFER METHODS 

3.1. Equation of Transfer and Definition of Coordinate 

System 

We follow Evans and Stephens (1991) and others in 
defining downward as positive, so that the direction of 
light propagation is defined by the vector {0, fi) or 
where 9 is the angle from the inward normal at the top 
of the atmosphere and yi = cos 9. The angle is the 
azimuth angle measured clockwise looking upward. We 
assume an unpolarized collimated incident solar beam 
in the direction (/ig, fio = 0). Light that is scattered 
backward toward the source then has the direction {—fJ-o-, 
tt). The Stokes vector I is defined as a 4-element column 
vector given by 


^total — fparaC^i^i^))'^ (1 fpara^Cr(S(l)) -\- 

^[cr(Q(0)) -F cr(Q(l))-H cr(Si(0))-H cr(Si(l))] (3) 

where we will hereafter refer to the last term as the Q 
cross section, as if it were due to a single transition, and 
its four individual contributions will be assigned the same 
wavenumber shift of 4161 cm^. 


2.3. Raman phase functions 

The generalized scalar phase function for scattering of 
unpolarized incident radiation by anisotropic molecules 
for all types of molecular scattering (including Raman 
and Rayleigh scattering) can be written as (Placzek 1959; 
Soris and Evans 1999): 


m 


3 r (1 -F d) -F (1 — <5) cos^ 9' 
21 2Td - 


( 4 ) 


I = 


I 

Q 

u 

V 


( 5 ) 


where I is the total intensity and (Q^ -\-U‘^ + V‘^)^ is the 
intensity of polarized light. Separate definitions for Q, 
U, and V are given by Hansen and Travis (1974). 

With these definitions, the vector radiation transfer 
equation can be written as 


dz 


“'SCcll.l. A 

Att j 


M(/x, 4>, p!, (f'fiiz, p.', (j)')dpL d(j)' -F 


'-1 


( 6 ) 


where 5 is the depolarization ratio, which is the ratio 
of perpendicular to parallel intensities scattered at 90°. 
The depolarization ratio is zero for isotropic particles 
(the pure Rayleigh case) and unity for isotropic scatter- 
ers. Eor pure rotational Raman scattering, 5 = 6/7 (Soris 
and Evans 1999), yielding P{9) = ^[13 -F cos^ 9], which 
has a peak deviation of only 7.7% from isotropy. De¬ 
polarization ratios for vibrational transitions for gases 
seem to be generally below about 0.38 (Bhagavantam 
1942). For the 4156 cm“^ Q-branch transition, mea¬ 
sured values range from 0.045 (Cabannes and Rousset 
1936) to 0.13 (Bhagavantam 1931). Either value implies 
a phase function that is closer to pure Rayleigh than to 
isotropic. However, the deviation from isotropic is great¬ 
est for the first scattering, for which the direct beam 
plays the largest role. Backscattered Q-branch light from 
the hrst scattering would be reduced by 30-33% using 
the assumption of isotropic scattering. But in sample 
calculations the reflected flux of Raman scattered pho¬ 
tons attributable to this transition is only 17% of the 
total (see Eig. [7]) and about half of that comes from in¬ 
cident light that already been diffusely scattered. Thus, 
the net effect of assuming isotropic scattering for the vi¬ 
brational transition will probably never exceed 15% of 
the Q-branch contribution to single scattering, and will 
generally be insignificant. Thus we will assume isotropic 
scattering for both rotational and vibrational transitions 
to benefit from the simplifications that it permits. 


where z denotes distance downward into the atmosphere, 
fcext and fcscatt are linear extinction and scattering coef¬ 
ficients, and where J is the source vector due to Raman 
scattering. Equation [6] is often written in terms of op¬ 
tical depth r, where dr = kextdz, and single-scattering 
albedo w = fcscatt/^ext- However, because Raman scat¬ 
tering transfers photons from one wavelength to another 
at a specific physical location in the atmosphere, the 
wavelength-independent distance ^ is a more convenient 
vertical coordinate than r. In Eq. [6] the scattering ma¬ 
trix M is a rotated form of the scattering phase matrix 
P, given by 

fi'= L(i 2 - 7r)P(cosa)L(ji) (7) 

where a is the scattering angle, ii is the angle between 
the scattering plane and the meridional plane of the in¬ 
coming ray, and 12 is the angle between the scattering 
plane and the meridional plane of the outgoing ray. Eor- 
mulas for the scattering and rotation angles, and the form 
of the rotation matrix, are given by Hansen and Travis 
(1974). In this formulation, the Stokes vector I includes 
direct beam as well as diffuse components of the radi¬ 
ation field, and thus there is no solar pseudo source as 
used by Evans and Stephens (1991). There is an implicit 
wavelength dependence in Eq. |6l but the wavelength de¬ 
pendence must be made explicit in the equation for the 








Accurate and Approximate Calculations of Raman Scattering. 


5 


source vector: 


ture points: 


1 MR 

J(z, n, ftnH 2 {z)am{vl) 

i=<d 

J Mr,(at, fl', ,z*,)dn', (8) 

where we use wavenumber iz = 1/X instead of wavelength 
because each Raman transition is associated with a fixed 
shift in wavenumber but a variable shift in wavelength. 
Here the quantity nH 2 (z) is the volume number density 
of H 2 molecules, and cross section for pho¬ 

tons at wavenumber = ly + Avi to excite transition i 
and exit at wavenumber iz, where hAvg/c is the energy 
change associated with the transition. The associated 
scattering matrix is Mr,. Note that the vector inten¬ 
sity l(z,/r', (/)', r'J') in the integrand is evaluated at the 
wavenumber of the incident photon, as are the cross sec¬ 
tion and scattering matrix in Eq. [5] If I is measured in 
photons, then the factor /, = 1. If I is measured in energy 
units, then /, = vlv*g. Although the radiation transfer 
equation we use will retain polarization in general, we 
will ultimately not include polarization in computation 
of Raman source contributions. 

There is also an extinction contribution from Raman 
scattering, which is given by 

riR UR 

kram = '^kraTn,i = '^nH2{z)aiu{v). (9) 

e=o e=o 

This is the sum over riR transitions of extinctions involv¬ 
ing photons of incident wavenumber i/ that experience 
Raman scattering and emerge as photons of wavenum¬ 
ber 1 / — Ai/,. The actual energy loss at wavenumber ly 
due to extinction is proportional to the incident energy at 
wavenumber v, but the source contribution at wavenum¬ 
ber V depends on the incident energy at other higher 
wavenumbers. 


3.2. Solution Method 

To solve the vector transfer equation we make use of 
Fortran vector radiation code developed and documented 
by Evans and Stephens (1991) and further validated 
by Sromovsky (2004) by comparisons with independent 
solutions of Sweigart (1970), Dlugach and Yanovitskij 
(1974), and Kattawar and Adams (1971). Using an 
order-m Fourier expansion of the intensity and phase ma¬ 
trix, Eq. iniis separated into 2m -I- 1 uncoupled equations 
that can be solved independently. The azimuthal compo¬ 
nents are separately solved and then combined to obtain 
intensity fields. This separation is valid for randomly 
oriented particles with a plane of symmetry, so that the 
16-element phase matrix has only six unique components 
(Hovenier 1969). We also assume that the incident radia¬ 
tion is symmetric in (j) and unpolarized. The zenith angle 
variation is discretized using a double-Gauss numerical 
quadrature. The radiance at any location in the atmo¬ 
sphere is thus represented as a vector involving three 
components: Stokes parameters, quadrature zenith an¬ 
gles, and azimuthal expansion mode. The intensity and 
source radiance vectors thus become vectors of ele¬ 
ments, where is the number of zenith angle quadra- 



Kdi) 


J(Mi) 

= 

Kd2) 

, J = 

J(M2) 






( 10 ) 


where each element or J(//i) is a vector of the form 
given in Eq. [S] The radiance field is separated into up¬ 
ward and downward hemispheres: U represents down¬ 
ward radiance (/i > 0) and l~ representing upward radi¬ 
ance (/r < 0). 

A model atmosphere is constructed by dividing the at¬ 
mosphere into uniform sublayers, each of which is char¬ 
acterized by an optical depth, a single scattering albedo, 
and scattering phase function. The code uses differential 
generators derived from Eq. |6]and makes use of doubling 
to compute reflection and transmission matrices for ho¬ 
mogeneous layers, which are then combined using stan¬ 
dard adding equations. The key equations are those of 
the interaction principle (Goody and Yung 1989): 


|+ = T+I+-(-R+Ii7+S+ 
|- =T-|--bR-|+-bS- 


( 11 ) 


where subscripts refer to the top (t) and bottom (6) 
boundaries of a layer, R and T are the reflection and 
transmission matrices, and S is the source vector for a 
layer, which equals Mz/fi for a differential layer. Using 
Eq. 111! it is straightforward to derive the properties of a 
combined layer (T + B) from the individual properties of 
the top (T) and bottom (B) layers (Evans and Stephens 
1991): 


■R 

■^T+B 

= R+ -H T+r+R 


fB 

Rt+b 

= R;^ -1- T^pT R 

B 

r+ 

T+ 

^T+B 

= TgT+T.^ 



T~ 

^T+B 

= T^r Tg 



C + 

^T+B ■ 

= s+-bT+r+(s+ 

-b 

R^! 

Sj+B ■ 

= Sj -i- T,pr (Sg 

-b 

Rb 

r+ = 

[1 — R^Rg] ^ 



r- = 

[1 — RgRiji] ^ 




( 12 ) 


Here the source adding equations refer to vector sources. 
The same equations can be shown to apply to the ma¬ 
trix sources that will be defined in a subsequent section. 
These can be converted to doubling equations by making 
the top and bottom properties equal. 

The first step in the solution is to convert Eq. [5] into 
a form that allows computation of the reflection and 
transmission matrices and the source vector. After az¬ 
imuthal expansion and zenith angle discretization, the 
vector transfer equation for azimuthal order m can be 
written in the form 

1+(z + Az) = T+ 1+(z) + R+1- (z + Az) + S+ , , 

I- (z) = T^l-(z + Az) + R-1+ (z) + S- ’ ^ 

where the matrices and source vector here refer to a dif¬ 
ferential layer of thickness Az. The single-scattering ap¬ 
proximation for the generator equations for reflection and 
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transmission matrices are then 




— (l 


flj 


) + 


kscatt^Z ^ |Mm(2:, V, ±^j>)\ii> 




IR^ 


'rnUjj'i' 


kscatt^Z |]V[^(z, Z/, zt/Tj , ) l^j/ (14) 

where i refers to Stokes vector index, j to zenith angle 
quadrature index, Wj' to the Gaussian quadrature weight 
value, and where Az is is usually chosen so that At = 
kextAz = 10“^, which Evans and Stephens (1991) have 
found to yield about 5 digits accuracy when using double 
precision. These expressions are different in form but 
equivalent to those given by Evans and Stephens. 

The differential generator for Raman source terms can 
be derived from the single-scattering form of Eq. |6] and 
discretization of Eq. [51 which results in the expression 
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1=0 

in which the individual matrix elements are given by 

fl X \krarn,^{z,vl)Az'^\MR^rn(.-^^J^J,±^JI■j')\ii' , (16) 

where the first ± superscript indicates the sign of and 
the second indicates the sign of /i', and where the Raman 
scattering coefficient per unit distance for transition £ is 
given by 

kramAz^^^i) = nHAzAmAl)- ( 17 ) 

Note that the source expressions depend on the 
wavenumbers of the incident photons {v^) that are 
shifted to the current wavenumber, while everything else 
in Eq. |T5] depends on the current wavenumber v. This 
is not a convenient form to work with because once the 
source vector is computed it applies to one specific inci¬ 
dent irradiance direction, meaning that a center-to-limb 
scan would require many repetitions of the solution algo¬ 
rithm. It is also unnecessarily complex because it retains 
the full generality of the phase function, which is not 
needed to accurately approximate the effect of Raman 
scattering. A simpler and more useful matrix formula¬ 
tion is developed in Section 13.31 
The solution is expressed in terms of reflection and 
source matrices for the entire atmosphere, from which 
intensities at desired irradiance and view angles are com¬ 
puted as described in Section 13.51 A uniformly spaced 
wavenumber grid is used so that the number of photons 
that are Raman scattered out of one bin will be trans¬ 
ferred to exactly one other bin. This requires that the 
wavenumber shift for each transition to be close to an in¬ 
tegral multiple of the grid spacing, as illustrated in Table 
4. Phase functions are introduced as Lagrange function 
expansions of an order that depends on the number of 
quadrature angles. Eor sharply peaked forward scatter¬ 
ing phase functions we avoid impractically large expan¬ 
sion orders, otherwise needed to avoid angular oscilla¬ 


tions in the reflected intensity, by employing the (5-Eit 
method of Hu et al. (2000). 

While the full polarization machinery of the Evans and 
Stephens (1991) code is retained in our modification, the 
Raman scattering, as we represent it, does not itself in¬ 
troduce any polarization. Further, because in most cases 
Rayleigh scattering is the dominant process in creating 
polarization effects, we usually make use of an approxi¬ 
mation that attributes all of the polarized intensity, given 
by (Q^ + U'^ + V^)^, to the Rayleigh scattering, and then 
use the algorithm described by Sromovsky (2004) to cor¬ 
rect the scalar calculation for polarization, rather than 
carry out the full vector calculation. The approxima¬ 
tion has been shown to be accurate to generally better 
than 1%, and entails no significant additional compu¬ 
tational burden, which otherwise would take about 40 
times longer. 

3.3. Matrix Formulation of the Source Function 

A matrix formulation of the source function is used to 
facilitate the computation of center-to-limb variations in 
reflectivity. Making the reasonable approximation that 
Raman scattering is isotropic (see Sec. 12.3L only the 
TO = 0 azimuthal component needs to include a Raman 
contribution and we can write 

fl X ^kraniAz^^'Zi)^z:^6aSi'i (18) 

and simplify Eq. [T5| to the form 

= J2'i=o^e0^iz,A) + l'(^)^l)) = 

E^=o^iGi'^A = 0,izf), (19) 

where the matrix Ge produces the sum of the upward 
and downward intensities at 2 by multiplication by the 
incident intensity at the top of the atmosphere (z = 0) 
at wavenumber lyf. There is an £ subscript on this ma¬ 
trix because it depends on transmission, reflection, and 
source matrices at the shifted Raman source wavenum¬ 
ber Ilf. Expressing the source in terms of the top-of- 
atmosphere incident irradiance vector allows us the com¬ 
putational convenience of expressing the source term as 
a matrix multiplier that has the same structure as the 
reflection matrix, which will be demonstrated in what 
follows. 

The Gi matrix can be derived from the interaction 
principle (Eq. Ill[) . The radiance at a given level in 
the atmosphere is proportional to the incident irradiance 
and can be expressed in terms of the properties of the en¬ 
tire atmosphere above that level (the top layer) and the 
properties below that level (the bottom layer). We as¬ 
sume here that the transmission out of the bottom layer 
is zero, i.e. if the atmosphere is thin, the bottom layer 
includes any surface reflection and/or absorption. The 
relevant source vector for the top layer Sj has multiple 
terms contributing from different Raman transitions that 
are each linearly related to the irradiance at any single 
wavelength (the source direction is the same at all wave¬ 
lengths). This can be shown by using the source addition 
equations given in Eq. 1121 If the sources for each layer 
to be combined are separately proportional to the top 
of atmosphere irradiance !(/ (= \'^{z = O^iif)), then the 
source for the combined layer will also be proportional 
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TABLE 4 

Wavenumber grid spacings that are commensurate with wavenumber shifts of primary Raman transitions. 


Transitions: 

Az/: 

S(0) 

354.39 cm-i 

5(1) 

587.07 cm-i 

Q + 5i 
4161.00 cm-l 



grid spacing Su 

Az//5z/ (integral steps) 

RMS error 

step sum 

12.240 cm-l 

28.953 (29) 

47.963 (48) 

339.951 (340) 

0.044 steps 

417 

13.640 cm-l 

25.982 (26) 

43.040 (43) 

305.059 (305) 

0.042 steps 

374 

17.780 cm-l 

19.932 (20) 

33.019 (33) 

234.027 (234) 

0.044 steps 

287 

39.260 cm-l 

9.027 (9) 

14.953 (15) 

105.986 (106) 

0.032 steps 

130 

58.620 cm-l 

6.046 (6) 

10.015 (10) 

70.983 (71) 

0.029 steps 

87 

118.860 cm-l 

2.982 (3) 

4.939 (5) 

35.008 (35) 

0.037 steps 

43 


NOTE: The step sum is proportional to the number of storage locations required to save informa¬ 
tion about shifted photons when calculating in sequence from shortest to longest wavelengths. 


to Iq". Thus we 
that satisfy 

s+K) = s+K)i+(i.;) Sb K) = s-(i.;)io+(z.;). 

( 20 ) 

Using these definitions, it is easy to show that 

= (1 -f Rb)(1 - R+Rb)-i [T+ -h S+R+Sb] + Sg, 

( 21 ) 

which defines the matrix at level ^ in terms of reflec¬ 
tion, transmission, and source matrices for the entire at¬ 
mosphere above z (subscript T for top layer) and that be¬ 
low z (subscript B for bottom layer), all at wavenumber 
. The source vector at z for the unshifted wavenumber 
V can then be written as 

s±(z.) = ^S,G,lo+(i.;) = ^ S,G,-^^lo+(z.) = Slo+(z.) 

( 22 ) 

in which the final two forms define a new matrix S and 
express the source contribution to the current wavenum¬ 
ber 1 / in terms of the solar irradiance at ly, using a scal¬ 
ing that is the ratio of solar spectral functions Fq{v) 
and Fq{vI), such that l+(r';) [Fq{vI)/Fq{v)\\+{ v). 

If the radiance field is measured in terms of energy per 
unit wavenumber, then Fq{v) is the solar energy spec¬ 
trum. If the radiance field is measured in photons per 
unit wavenumber, then Fq{v) is the solar photon spec¬ 
trum. Substituting the expression given by Eq. 1 181 yields 
the computationally useful form 


tational procedure described below. 

3.4. Computation of the Source Matrix 

Atmospheric layers that have uniform extinction co¬ 
efficients and phase functions become nonuniform due 
to the vertical variation of the Raman source contribu¬ 
tion. This requires us to use a much larger number of 
layers than would be necessary were Raman scattering 
not considered. We typically use 30-80 logarithmically 
spaced pressure levels between 0.0003 and 100 bars, with 
80 levels required for accuracy of a few tenths of 1%. 
The accuracy achieved in any particular problem needs 
to be assessed by running test calculations using a larger 
number of levels. 

Each Raman solution requires a spectral calculation. 
Following Betremieux and Yelle (1999), we start at the 
highest wavenumber and work downward so that mul¬ 
tiple Raman scattering is automatically included in a 
single pass. At the highest wavenumber, the radiation 
transfer problem is solved with no source contribution, 
but with Raman extinction included. At each wavenum¬ 
ber including the first, the G matrix is computed at each 
layer bottom boundary L, using Eq. The algorithm 
then computes for each transition each layer L, and 
each possible incident quadrature index j", the function 

^ I. - 

(24) 

3' 


can define source matrices and Sg 


\- ^ A V 0 


j' 


J l^XjFQ{v)i 


(23) 


in which we have here made a specific choice that FQ(y) 
is a solar spectrum in energy per unit wavenumber, and 
where i" = 1 is the only relevant value of the third in¬ 
dex for the assumed unpolarized incident irradiance, and 
f = 1 is the only relevant value of the first index for 
the assumed isotropic and Raman scattering phase func¬ 
tion, according to Eq. [TH The ± superscript on S was 
dropped, and never used on S, because both values are 
the same for a differential layer (see Eg. IT^ . Note that 
the first bracketed factor depends on vf while the second 
bracketed factor depends only on wavenumber v and z, 
and is independent of £. This form suggests the compu- 


which is essentially the same as the left bracketed factor 
in Eq. [23] divided by Az. This is essentially the photon 
loss term per unit distance for each layer and each tran¬ 
sition. Here everything is evaluated at the same wave¬ 
length and thus the G matrix needs no i dependence. 
This form uses the average of G at the top and bottom 
of each layer to represent the radiation field everywhere 
within the layer. 

As the computation proceeds to longer wavelengths, 
the source contributions at the current wavenumber are 
then obtained by reading the stored values of F[ and 
computing the source matrix as follows: 


(S^)ijij'' 


^ i 
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Az. 


(25) 


For each wavenumber at which values of FI are read, 
the values for transition i will have been stored during 
calculations made at wavenumber n + Ani. 
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Once H is read from storage and the matrix is 
computed, the computation of the radiation field pro¬ 
ceeds with the assumption that the source is uniformly 
distributed over the layer L. This is necessary to make 
use of the doubling equations for a uniform source, but is 
only a good approximation for layers that are relatively 
optically thin at the wavenumber of the incident photons. 
While the input source is assumed to be uniform through¬ 
out each layer, and is input as the differential source after 
dividing by the physical thickness of each layer, the dou¬ 
bling process does modify the source function in accord 
with the absorption and scattering processes that take 
place at the wavelength of the scattered Raman photon. 


3.5. Computation of Intensity 

After computing the source matrix for the entire at¬ 
mosphere using the same adding equation that is given 
for the source vector in Eq. we then compute the 
outward intensity at the top of the atmosphere. For ar¬ 
bitrary directions of incidence and reflection this is for¬ 
mally expressed as 


lo(Ai»= R- + S (/r,^o,</'-0o) 


I^qFo/tt 


IJ-qFo/t^ 

0 

= R- 

0 

0 

0 

0 


0 


(26) 


where the incident flux is Fq in direction (0o, 4>o) and the 
combined matrix S“ uses the minus superscript because 
direction does make a difference for the inhomogeneous 
layer consisting of the entire atmosphere. By converting 
the source function into a matrix we have thus put the 
Raman contributions in the same form as the reflection 
contributions, allowing us to compute a pseudo reflection 
matrix R~ and use the same mathematical machinery to 
compute spatially resolved intensity profiles. 

The azimuthal expansion of R“ is given by 

= Rq (^,/To) 
oo 

{fi, fio) cos{m{(l) - (fo)) (27) 

m—1 


Where sine components of the expansion are zero by the 
assumed symmetry of the incident radiation field. Since 
our solution only provides R“ at quadrature points, we 
carry out a Legendre polynomial interpolation to obtain 
intensities at other angles. 


3.6. Validation 

How do we know that our computations are correct? 
The basic radiation transfer algorithm without Raman 
scattering has been compared with independent calcu¬ 
lations for selected cases and shown to be in excellent 
agreement. This was done by Evans and Stevens (1991) 
for their original code and by Sromovsky (2004) for the 
modification used here (as noted in Sec. 13.2p . but nei¬ 
ther of these sets of comparisons deals with the Raman 
scattering additions, which are harder to validate. A 
few comparisons with past calculations can be made, but 
most of these lack the resolution, relevance, or rigor that 
are needed to serve as appropriate standards of accuracy. 


They do provide a reasonable sanity check, however, and 
some useful comparisons are made in Section 15.41 An¬ 
other validation is to show by means of test cases that the 
algorithm satisfies conservation of photons. This is done 
for the case of a monochromatic incident flux in Section 
15.11 Conservation is unlikely to be achieved if the Ra¬ 
man transfers are improperly computed. A last valida¬ 
tion, though an indirect one, is to demonstrate that the 
algorithm can reproduce features in the observed sp ectr a 
with reasonable atmospheric structure models ISec. 15.51) . 

4. ATMOSPHERIC AND SOLAR MODELS 
4.1. Model Atmosphere of Neptune 

Neptune’s thermal structure is obtained from Voyager 
radio occultation observations at ingress (45° S) and 
egress at 61° N). We used the Hinson and Magalhaes 
(1993) analysis for p < 1 bar and the Lindal (1992) re¬ 
sults for 1 bar < p < Q bars, with an offset of 1.0 K 
added to match the the Hinson and Magalhaes profile at 
1 bar. For p > 6 bars, we extrapolated using the nearly 
adiabatic lapse rate at the 6 -bar level (« -0.94 K/km). 
The two temperature profiles differ insignificantly except 
in the stratosphere. The lower egress temperatures in 
the tropopause region reduce the CH 4 mixing ratio and 
opacity to a small degree that is only noticed near 0.89 
/rm, where a 10 % increase in geometric albedo is pro¬ 
duced. We assume a fixed profile at all latitudes for ra¬ 
diation transfer analysis. Our altitude scale is computed 
using gravitational acceleration at 45° S {g = 11.1 m/s^). 
The Hinson and Magalhaes profile is derived assuming a 
gas composition of 81% H 2 and 19% He (Conrath et al. 
1991). It is possible that N 2 may also be present at a 
mixing ratio as high as 0.3% (Conrath et al. 1993), in 
which case our altitude scale would be modified some¬ 
what. The Lindal profile assumed 2% CH 4 in the tropo¬ 
sphere, which is close to the currently accepted value of 
2.2% (Baines et al. 1995) that we use for computing CH 4 
opacity. The tropospheric mixing ratio is consistent with 
CH 4 condensation at 1.4 bars, above which we assume 
the saturated vapor pressure until we reach the strato¬ 
sphere, at which point the mixing ratio is the smaller of 
the saturated mixing ratio and the stratospheric limit of 
3.5 xl0“'* (Baines and Hammel 1994). 

The pressure at which the one-way vertical optical 
depth reaches unity for each opacity contributor, and for 
the total of all contributors is illustrated in the top panel 
of Fig. m The levels at which the total optical depth 
reaches values from 0.3 to 100 are displayed in the bot¬ 
tom panel. The penetration depth of sunlight into Nep¬ 
tune’s atmosphere is limited by Rayleigh scattering at 
short wavelengths and by CH 4 and H 2 collision-induced 
absorption (CIA) at long wavelengths. The deepest pen¬ 
etration is at 0.935 /rm, where there is a relatively clean 
CH 4 window, a window in the CIA spectrum, and a low 
Rayleigh cross section. We used CH 4 absorption coeffi¬ 
cients derived by Karkoschka (1994) from planetary ge¬ 
ometric albedo observations using a technique described 
by Karkoschka and Tomasko (1992). Karkoschka esti¬ 
mates a 5% uncertainty in his 1994 absorption coeffi¬ 
cients, plus an additional uncertainty in the continuum 
baseline: 0.0003 km-am“^ at 400 nm and 0.02 km-am“^ 
at 1000 nm (a factor of 2 every 100 nm). This con¬ 
tinuum uncertainty leads to uncertainties in the CH 4 
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Fig. 3.— Vertical temperature structure of Neptune based on 
ingress (solid) and egress (triple dot dash) occultation observations, 
and volume mixing ratios of H 2 (dashed), He (dotted), and CH 4 
(dot dash). The CH 4 mixing ratio is 3.5 xl0“^ in the stratosphere, 
0.022 in the deep troposphere, and otherwise follows the saturation 
vapor pressure curve. The horizontal dotted line at 1.43 bars in¬ 
dicates the point at which methane reaches the saturation vapor 
pressure in the troposphere. 


window regions that are important to determination of 
the scattering properties and pressure levels of clouds 
in the several bar range. There is also a likely bias er¬ 
ror for weak CH 4 bands superimposed on Raman scat¬ 
tering effects due to the nonlinear relationship between 
I/F and single-scattering albedo (this is discussed in Sec. 
E3]). The CIA values are obtained by interpolating ta¬ 
bles of pressure and temperature dependencies provided 
by Alexandra Borysow, and available at the Atmospheres 
Node of NASA’S Planetary Data System. The average 
Rayleigh scattering cross section per molecule was com¬ 
puted using the equation 


CTRay 


3A47V2 “ ) 6-7Si 


(28) 


from Hansen and Travis (1974), where Vi is the volume 
mixing ratio of the ith gas, Ug^i is its refractive index, 
and 5i is its depolarization factor. We used depolariza¬ 
tion ratios of 0.0221 for H 2 and 0.025 for He (Penndorf 
1957; Parthasarathy 1951) and refractive index values 
from Allen (1964). Because no depolarization values were 
give for CII 4 , we used CO 2 values. We also used the 
wavelength dependence of ammonia’s refractive index to 
approximate that of CH 4 . The error in these latter ap¬ 
proximations is not significant because of the small CH 4 
mixing ratio. 


4.2. Cloud and Haze Models 

Two model structures with aerosols were chosen to il¬ 
lustrate the effects of stratospheric haze and lower alti¬ 


Fig. 4.— Top: Pressure level at which Rayleigh scattering 
(dashed), Raman scattering (dotted), CH 4 absorption (triple dot- 
dash), H 2 CIA (dot-dash), and the sum of all opacities (solid) 
accumulate a one-way vertical optical depth of unity, plotted ver¬ 
sus wavelength. Bottom: Pressure at which one-way total vertical 
optical depth reaches 0.3, 1, 3, 10, 30, and 100, plotted versus 
wavelength. 

tude cloud aerosols on Raman scattering features and to 
test approximations under more realistic conditions than 
provided by a clear atmosphere. Both models contain a 
high-altitude absorbing haze of Mie-scattering spherical 
particles with a Hansen (1971) gamma size distribution 
of effective radius a = 0.2 /im and variance b = 0 . 02 , 
where the relative number per unit radius interval at ra¬ 
dius r is proportional to r*^4-3f))/b exp(—r/ [ah)). The par¬ 
ticle size is from Pryor et al. (1992). The refractive index 
of the haze is assumed to have a real value of 1.44 and 
A-dependent imaginary values given in Fig. [SJ which dif¬ 
fer from the variation inferred by Courtin (1999). Other 
index variations could also have been used, with com¬ 
pensating changes in other aerosol components. Finding 
a tightly constrained fit is left for future work. While 
Courtin’s r = 0.1 haze extends downward only to 20 
mb, our haze has twice the optical depth and extends 
much more deeply; it has a uniform mixing ratio between 
100 and 800 mb, which is more similar to the distribu¬ 
tion inferred by Moses et al. (1995). Both of our haze 
models also contain a 3.8-bar cloud of isotropic scatter- 
ers, with a single-scattering albedo 0.99. This cloud is 
at the level that is often considered to be the top of a 
semi-infinite cloud (Baines et al. 1995). Haze Model I 
has this deep cloud set to optical depth 0.5. For Haze 
Model H this cloud is set to unit optical depth and a 
second cloud is placed at 1.3 bars, the approximate level 
expected for the base of a CH 4 ice cloud assuming a 
CH 4 mixing ratio of 2.2%. For this cloud we assumed 
r = 2, a single-scattering albedo of 0.99, and a double 
Henyey-Greenstein phase function of the form P{0) = 
fiPiid) + (1 - fi)P2i0), with 51 = 0.9, g 2 = - 0 . 11 , and 
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/i = 0.42, where P,( 6 ') = {I-gf) / {l+gf-2gicos{0)f/'^. 
The phase function parameters are due to Pryor et al. 
(1992), based on high-phase angle Voyager images. Our 
haze and cloud structures are not meant to be optimum 
fits to Neptune’s geometric albedo spectrum, but rather 
a sampling of possible structures that might be encoun¬ 
tered. Nevertheless, the part of the spectrum below 0.5 
/rm is a fairly good match to the obse rved geometric 
albedo spectrum, as demonstrated in Sec. 15.51 The Haze 
11 model mainly serves as a test case for approximations. 



Fig. 5.— Upper haze imaginary index assumed here (solid curve) 
compared to that inferred by Courtin (1999) (dashed curves) and 
by Pryor et al, (1992) (plotted point). 


4.3. Models of the Solar Irradiance Spectrum 

A key input to modeling of Raman scattering is 
an accurate solar spectrum of the appropriate spec¬ 
tral resolution. For A > 0.41 pm we used a model 
spectrum (Kurucz 1993) normalized to the 0.41 to 
0.87 pm results of Neckel and Labs (1984) and for 
0.12 pm < A < 0.41 pm we used measurements by 
the Upper Atmospheric Research Satellite using the 
SOLSTICE instrument (Woods et al. 1993). The 
UARS spectrum was obtained from the UARS web site 
(ftp://rescha.colorado.edu/pub/solstice/sol_hires_ 200 .dat) 
It has a nominal resolution of 0.2 nm, but is actually 
closer to a resolution of 0.5 nm, based on comparisons 
with convolutions of the very high resolution Kurucz 
model spectrum. We created our standard solar ref¬ 
erence spectrum (Fig. I®]) by convolving the combined 
spectrum with a triangular sampling function to obtain 
a nominal FWHM resolution of 0.35 cm“^ sampled at 
0.1725 cm“^. The nominal wavelength resolution varies 
from 0.14 nm at 0.2 pm to 3.5 nm at 1 pm (the resolution 
is 1 nm at 534 nm). Although this has nominally a 
uniform wavenumber resolution, it is actually limited 
by the UARS observations to about 0.5 nm at the 
shortest wavelength, where it is oversampled by a factor 
of 2 or more. The spectrum doesn’t actually reach the 
nominal resolution until the transition to the Kurucz 
spectrum at 0.41 pm. The solar spectrum matches the 
Karkoschka ground based observing resolution of 1 nm 


at 0.534 pm, but is 3.5 times worse near 1 pm. Thus 
careful comparisons between observation and model 
calculations need to adapt to differing resolutions. 

Following Courtin (1999), we selected a UARS spec¬ 
trum from 29 March 1992 to provide optimum compat¬ 
ibility with the FOS observations from 19 August 1992. 
Courtin found a peak-to-peak difference of only 2% be¬ 
tween two low-resolution UARS daily spectra obtained 
on 2 June 1992 and 19 August 1992. Thus, solar vari¬ 
ability is unlikely to be a major error source in com¬ 
paring calculations with FOS observations. Most of the 
solar variability is restricted to A < 0.26 pm, so that 
comparison with the groundbased observations made in 
1993 (Karkoschka, 1994) are even less influenced by solar 
variability. 

5. CHARACTERISTICS OF RAMAN SPECTRA 

5.1. Raman Scattering for a Monochromatic Source 

The behavior of monochromatic incident photons pro¬ 
vides useful insight into the workings of Raman scat¬ 
tering in Neptune’s atmosphere and a useful validation 
of our basic computational algorithm. For this exam¬ 
ple we use the three Raman transitions given in Table 
4, a wavenumber spacing of 58.62 cm“^, an atmosphere 
with CH 4 absorption and H 2 CIA, an equilibrium dis¬ 
tribution of H 2 , but no aerosols. The propagation of 
photons from the injection wavelength of 0.228247 pm 
(43812.2 cm“^) is illustrated in the bottom panel of Fig. 
Ul The vertical scale in this plot is the ratio of the total 
reflected flux at normal incidence to the incident photon 
flux. The first peak shows that 54.6% of the incident 
photons are reflected by the atmosphere at the incident 
wavelength. The remaining 45.4% of the photons are 
shifted to longer wavelengths by Raman scattering. The 
second peak, containing about 8 . 1 % of the photons is due 
to the 5'(0) transition. This second peak is the fraction of 
single-scattered 5'(0) photons that exit the atmosphere; 
but many S(0) photons are further shifted to even longer 
wavelengths. The third peak, containing 2.2% of the in¬ 
cident photons, is due to the <S'(1) transition. Many addi¬ 
tional peaks are due to multiple Raman scattering with 
various combinations of S'(l) and S{0). 

The relatively large peak near 0.252 pm, containing 
8.1% of the incident photons, is the first due to the Q 
transition. This is the same size as the 5'(0) peak and 
thus seems surprisingly large given that the Q transition 
has about half cross section of the 5'(0) transition (the ra¬ 
tio is 0.541 at 0.228 /xm). However, the 5'(0) contribution 
is multiplied by fpara, which is below 0.75 in the pres¬ 
sure range likely to contribute most at this wavelength 
(Fig. [T]). One might thus expect the Q contribution to 
be at least 72%, of the 5'(0) contribution on this basis 
alone. However, there are other factors also at work, 
such as multiple scattering and the vertical distribution 
of Raman photons. In fact, if the Rayleigh scattering 
cross section is made wavelength-independent, then the 
Q peak becomes only about 57% of the S'(O) peak, which 
is less than expected from cross section and fpara consid¬ 
erations. Apparently, in the real atmosphere the decreas¬ 
ing opacity of the atmosphere with increasing wavelength 
allows more of the Raman photons to make it out of the 
atmosphere than would otherwise be the case. Because 
the Q photons undergo the largest wavelength shift, they 
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Fig. 6 .— Solar irradiance spectrum in energy flux units per wavenumber, using the UARS spectrum (Woods et al. 1996) for wavelengths 
less than 0.41 /rm, and a normalized model spectrum for longer wavelengths (Kurucz, 1993b), convolved to a wavenumber resolution of 
35 cm~^ and sampled at twice that frequency. The letters C-K label the more prominent Fraunhofer absorption lines. At the bottom are 
scaled and offset ratio plots illustrating fill-in features (smoothed spectrum/unsmoothed), 5(0) ghosts (unsmoothed/smoothed with 354.4 
cm“^ wavenumber shift), 5(1) and Q ghosts (unsmoothed /smoothed with 587cm~^ and 4161 cm~^ shifts respectively). 


exhibit the largest effect. 

The conservation of photons illustrated in the top 
panel of Fig. [7] is computed as follows. For an incident 
unpolarized flux Fq{v) at zenith cosine /iw, the quadra¬ 
ture value closest to unity, the total reflected flux at v is 
given by by the discrete summation 


Fr{v) = (v) 


klNl 


(29) 


where, for a monochromatic incident flux at iz = izo, only 
the reflection term contributes v = vq and only the 
source term contributes at iz < izp. Because of the ma¬ 
trix formulation of the source term that is used here, we 
must introduce a pseudo-monochromatic flux that has a 
negligibly small constant offset for iz izq, to avoid divid¬ 
ing by zero in Eg. 1251 Using discretization of wavelength 
as well, in which the incident wavenumber is izq and 
Vj = izg — j Azz, we can write the pseudo-monochromatic 


incident flux as F^fivj) x (5o,j + e) implying a photon 
line flux (flux at wavenumber vq) of Fq{vq)/( hv^). The 
cumulative fractional photon flux can then be written as 

) — Fo(i/o) ~ ^ Vj ~ 

(30) 

where J < v^jIS.v and where we ignored e compared to 
unity in the first term and made use of the fact that 
the source term is zero at v = jzq. Although the sec¬ 
ond term has an explicit multiplication by e, it is not 
negligibly small because the source term itself contains 
a division by e. For a conservative atmosphere that is 
either semi-infinite or bounded by a unit-albedo surface, 
/cum should evaluate to unity if the sum is extended to 
the last wavenumber. The cumulative flux shown in Fig. 
[7] falls somewhat short of unity because of CH4 absorp¬ 
tion and H 2 CIA. 
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Fig. 7.— Bottom: Photon flux ratio vs wavenumber for an incident monochromatic flux at 43812.2 cm“^ (228.247 nm). Selected flux 
peaks are labeled by Raman transition or combination of transitions that produced them. Middle: extended version of bottom panel. Top: 
Cumulative fractional photon flux. 


The middle panel of Fig. [7] shows that Raman photons 
can undergo a surprisingly large number of scatterings 
before they leak out at the top of the atmosphere. In 
fact, as shown by the plot of cumulative flux in the up¬ 
per panel, even by 0.45 /xm, 1.6% of the incident photons 
are still inside the atmosphere. This occurs because the 
atmosphere is conservative at short wavelengths and be¬ 
cause photons scattered deep within the atmosphere are 
more easily lost to another Raman scattering than lost 
to transmission out of the atmosphere (to space). Af¬ 
ter each additional Raman scattering, the optical depth 


to space gets smaller, as does the cross section for fur¬ 
ther Raman scattering (Fig.[2), so that photons are more 
easily lost to space at longer wavelengths. 

The vertical distribution of Raman scattered photons 
is illustrated in Fig. [8] for the Q transition, which pro¬ 
duces peaks spaced by 4161 cm“^ (Fig. [7]). The curve of 
highest magnitude is due to one scattering, the second 
due to two scatterings, and the last is due to eight suc¬ 
cessive scatterings involving the Q transition. Because of 
the large wavenumber shift of this transition, these peaks 
decline more slowly than those related to the S'(O) and 
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S'(l) transitions. The peak at the shortest wavelength 
is at a location where the incident radiation field times 
atmospheric density is at a local maximum, which occurs 
near 300 mb. After each additional scattering, the Ra¬ 
man photons diffuse both upward, where they leak out 
of the atmosphere, and downward, where they provide a 
reservoir for later contributions. It is worth noting that 
some of the Raman scattered photons that are trans¬ 
ferred to a longer wavelength can leave the atmosphere 
directly without any further scattering (those at high al¬ 
titude and directed upward). The result is a movement 
of the peak source intensity further downward after each 
successive scattering. 



Source/km (j= 4, j"= 4) 


Fig. 8. — Relative source function intensity per km ver¬ 

sus geometric mean pressure of the layer boundaries, for selected 
wavelengths (labeled by nanometers), for a clear Neptune atmo¬ 
sphere. 

Because of vertical variation in the source function, it 
was necessary to use 72 logarithmically spaced gas layers 
to avoid cumulative photon fluxes exceeding unity. When 
cross sections were made wavelength independent, but 
evaluated at 0.4 ^m, the accuracy for a given number 
of levels was much higher. Monochromatic calculations 
were run with three different choices for the number of 
logarithmically spaced atmospheric levels between 0.0003 
and 100 bars: 18, 36, and 72. The ratios to the 72-level 
fluxes, were 0.991-0.998 for the 18-level result, and 0.998- 
0.9996 for the 36-level result. Thus it seems likely that 
errors for 72 levels probably do not exceed 0.1%, and 36 
levels are probably within 0.2%. The ratios to the 72- 
level cumulative flux at 1 /im were 0.9976 for 18-levels 
and 0.9995 for 36 levels. 

5.2. Raman Scattering of a Flat Spectrum with an 
Absorption Line 

For an irradiance spectrum with a A-independent pho¬ 
ton flux, one might expect the apparent geometric albedo 
with Raman scattering to be the same as without it. 
One could argue that for every photon lost from a given 
wavenumber, there is one gained from a higher wavenum¬ 


ber. This argument works well if atmospheric properties 
are also independent of wavelength. But in a real at¬ 
mosphere, Raman and Rayleigh scattering cross sections 
vary with wavelength, so that the vertical distribution 
and fraction of photons lost will be somewhat different 
from the vertical distribution and fraction of photons 
gained, even with a A-independent incident photon flux. 
In spite of this complication, the results shown in Fig. 
are roughly consistent with this argument. The calcula¬ 
tion shown here is for a clear Neptune atmosphere with 
CH 4 absorption and CIA turned off, and exposed to an 
irradiance spectrum that has a constant photon flux ex¬ 
cept for a pseudo-solar absorption feature at 0.311641 
fj,m (32088 cm“^), where the photon flux was dropped 
to 20% of its value elsewere. The irradiance plotted here 
is the spectrum of energy per unit wavenumber, normal¬ 
ized to unity at 0.2 /rm, which corresponds to a flat pho¬ 
ton spectrum. Note the considerable distance from the 
starting wavelength to the wavelength at which the ge¬ 
ometric albedo approaches the Raman-free semi-infinite 
H 2 Rayleigh value of 0.7908 (Sromovsky 2004). That 
is not surprising because of the importance of multiple 
Raman scattering evident from Fig. 0 Calculation were 
run with surface albedos of 0 and 1 , to show that for 
A > 0.6 ^m even 100 bars of atmosphere is not enough 
to obscure the lower boundary when only Rayleigh and 
Raman scattering are considered. With CH 4 absorption 
turned on, the effects of the bottom boundary are no 
longer apparent. 

Among the sharp spectral features present in the re¬ 
flection spectrum, the largest peak is at exactly the wave¬ 
length of the absorption spike in the irradiance spec¬ 
trum. This peak arises because in addition to photons 
reflected at the same wavelength, there are also photons 
Raman shifted from shorter wavelengths. Because more 
are shifted in than shifted out, sharp absorption features 
like this get partially filled in the reflected flux. When the 
reflected flux is then divided by the irradiance to obtain 
the albedo, the feature stands out as a positive spike. So¬ 
lar absorption features like this are the most prominent 
features in the Raman spectra of the outer planets. The 
other, much smaller downward spikes are ghosts of the 
irradiance spectrum, displaced by the wavenumber shift 
of the Raman transition. In the detailed view shown in 
the lower part of Fig.|9]the ghost features are labeled by 
the transitions that produced them. Note the relatively 
small size of the ghost features. The S'(O) ghost feature 
has an amplitude of about 5.8% of the geometric albedo, 
while the corresponding absorption feature in the irradi¬ 
ance spectrum has an amplitude of 80%; the S'(l) and Q 
features have amplitudes of about 2.7% and 4.8% respec¬ 
tively. The near equality of Q to S{0) ghost amplitudes 
in spite of their very different cross sections has the same 
explanation as the relatively large Q contribution peak in 
the monochromatic spectrum, discussed in the previous 
section. 

5.3. Raman Scattering of the Solar Spectrum 

The effect of Raman scattering on Neptune’s aerosol- 
free geometric albedo is illustrated in Fig.fTOl where spec¬ 
tra including polarization are shown for three calcula¬ 
tions of increasing complexity, first with only Rayleigh 
scattering, then including also molecular absorption, and 
finally including Raman scattering. The calculations 
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Fig. 9.— Upper: Flat photon irradiance spectrum (dot-dash) 
with an absorption feature at 0.311641 fim (32088 cm“^) and cor¬ 
responding geometric albedo spectra of Neptune for unit (dashed) 
and zero (solid) surface albedos. The dotted line indicates the geo¬ 
metric albedo for a Raman-free semi-infinite Rayleigh atmosphere 
(p=0.7908). Lower: Geometric albedo vs wavenumber, showing 
Raman transition assignments for the main ghost features. 

were started at 0.2 ^m, and are accurate to better than 
1% at wavelengths beyond 0.22 /im. Experiments with 
different starting wavelengths show that an accuracy of 
1 % or better can be achieved at a given short wavelength 
by starting the calculation at a corresponding wavenum¬ 
ber that is higher by twice the wavenumber shift of the 
Q-branch transition. This relatively small starting off¬ 
set is effective because of the steep gradient in the solar 
spectrum at short wavelengths (see Fig. IH]). The calcula¬ 
tion was carried out with a step size of 17.78 cm“^, using 
a solar spectrum with a nominal resolution of 35 cm“^. 
With only conservative Rayleigh scattering, Neptune’s 
geometric albedo would be a wavelength-independent 
0.791 (Sromovsky 2004). Adding methane absorption 
and collision-induced hydrogen absorption (CIA) only re¬ 
duces the albedo at wavelengths beyond 0.4 ^m. Raman 
scattering produces profound effects throughout the UV 
and most of the visible spectrum. 

Besides the introduction of sharp spectral features, Ra¬ 
man scattering also reduce Neptune’s baseline geometric 
albedo by nearly 25% in a clear atmosphere. In Fig. (TU] 
the albedo decreases (shown as red) appear mostly at 
shorter wavelengths, but also appear in the vicinity of 
most of the local peaks in the near-lR reflectivity. The 
albedo gains (shown as black) occur not just at the deep 
minima in the solar spectrum, but also at the absorp¬ 
tion maxima in the CH 4 spectrum, which correspond to 
minima in the reflectivity spectrum. 

We find that the fill-in of the reflectivity minima in 
the near IR is only about 4% of the geometric albedo. 
This is in sharp contrast to the Cochran and Trafton 
(1978) conclusion that in the cores of the strong CH 4 
bands at 0.86 /xm, 0.89 fim and 1.0 /xm, nearly all of 


the photons that leave the atmosphere have been Ra¬ 
man scattered. However, their conclusion is a result of 
what we now know to be inappropriate assumptions con¬ 
cerning the amount and vertical distribution of methane. 
Their standard model had a constant mixing ratio of 
0.005 relative to H 2 throughout the atmosphere. This 
resulted in CH 4 supersaturation for p < 1 bar, by fac¬ 
tors of ^100 at the tropopause and more than 10 in the 
stratosphere. This resulted in a very low single scatter¬ 
ing albedo of a; ~ 0.002 at 0.89 /xm, which even for a 
semi-infinite atmosphere would yield a very low geomet¬ 
ric albedo of only 0.00036, using the single-scattering ap¬ 
proximation p{uj) = 0.1855 X w for w <C 1 (Sromovsky 
2004). That is about 20 times smaller than the geometric 
albedo we computed for our assumed CH 4 proHle, which 
never exceeds saturation levels and is limited to 3.5 x 10“"'’ 
in the stratosphere. Our distribution leads to much more 
scattering at the center of the disk and significant limb 
brightening (Sromovsky 2004), while the uniform mixing 
ratio of Cochran and Trafton makes the atmosphere dark 
at all view angles. 




Fig. 10.— Upper: Neptune’s geometric albedo spectrum for mod¬ 
els without aerosols assuming conservative Rayleigh scattering with 
polarization (green), assuming Rayleigh plus molecular absorption 
(red), and assuming Rayleigh scattering, molecular absorption, and 
Raman scattering (black). Lower: Spectral differences due to Ra¬ 
man scattering, with Raman losses (red) and Raman gains (black) 
compared to 0.04 X the geometric albedo spectrum (green). 

5.4. Comparison with other Raman Calculations 

The earliest accurate calculation of a Raman spectrum 
relevant to Neptune was by Cochran and Trafton (1978) 
for a pure H 2 atmosphere assuming an equilibrium pop¬ 
ulation of ortho and para states. Their most relevant 
case is for 500 km-amagats of H 2 over an opaque cloud 
layer (of unspecified but apparently high albedo). That 
amount of H 2 corresponds to an atmospheric depth of 
approximately 6 bars, which is deep enough to approach 
semi-infinite behavior only for A < 0.34 ^m. At 0.5 /xm 
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the surface albedo can make a difference of 0.25 in geo¬ 
metric albedo. Using 2/3 of the Ford and Browne (1973) 
Raman cross sections to approximate the values used by 
Cochran and Trafton, we computed geometric albedo val¬ 
ues for a He-free atmosphere with a unit-albedo Lam¬ 
bertian surface placed at the 6-bar level, and ignored 
polarization. This spectrum, smoothed to a resolution 
of 5-nm, is displayed in Fig. [TT] where it compares well 
with the even lower resolution spectrum of Cochran and 
Trafton. 

The other prior calculation displayed in Fig. [TT]is due 
to Savage et al. (1980), who used the same iterative 
computational scheme first presented by Cochran and 
Trafton (1978). According to Savage et al. , they used 
the Ford and Browne (1973) Raman cross sections that 
we used, but ignored polarization and assumed a pure H 2 
atmosphere, which was deep but of unspecified depth. 
In spite of their use of different Raman cross sections, 
the Savage et al. and Cochran and Trafton results agree 
with each other and with our calculations using 2/3 of 
of the Ford and Browne cross sections, raising questions 
about what cross sections were actually used by Savage 
et al. (1980). Our calculation with the Ford and Browne 
cross sections, and other conditions being the same, are 
shown as the thinner dot-dash curve in Fig. [TTJ Clearly, 
this spectrum is not compatible with the Savage et al. re¬ 
sults. All these calculations have a stronger upward slope 
than we obtained for a truly deep atmosphere (Fig. llT))) 
because the surface reflection becomes more visible at 
longer wavelengths. This can be seen from the difference 
between calculations for zero and unit surface albedos, 
which are shown in the figure. Other effects occur when 
we carry out more realistic calculations. The presence 
of helium dilutes the Raman absorption somewhat, pro¬ 
ducing about a 2% increase in the baseline value; us¬ 
ing the correct Raman cross sections causes an albedo 
decrease of about 0.02-0.03 (a 3-5% decrease); and po¬ 
larization increases the geometric albedo by about 0.04 
(Sromovsky 2004) (a 6.7% increase). The approximate 
net result is a very small difference between the accurate 
calculation with polarization and the earlier H 2 -only cal¬ 
culations without polarization. 

5.5. Comparison with Observed Spectra 

The International Ultraviolet Explorer (lUE) made 
full-disk observations of Neptune during 2-5 October 
1985 (days 275-277) at a phase angle of 1.87°and a sub¬ 
earth latitude of-19.95°planetographic (-19.325°centric). 
Wagener et al. (1986) converted these observations to ge¬ 
ometric albedo assuming an equatorial radius of 25240 
km and an oblateness of 0.021, which are average stellar 
occultation results referring to a pressure of 1 microbar 
(French 1984). Using measurements of 24764 km and 
0.01708 (Davies et al. 1992), which refer to the 1-bar 
level, we obtain a conversion factor of 1.035 for adjust¬ 
ing the Wagener et al. albedo values (this disagrees with 
the factor of 1.047 used by Courtin (1999)). We did not 
make a correction for phase angle, though that might in¬ 
crease their values by as much as 1% if the phase curve 
at UV wavelengths is similar to that observed at longer 
wavelengths (Sromovsky et al. 2003). The adjusted lUE 
results are compared with our clear-sky and Haze I Nep¬ 
tune spectra in Fig. where error bars indicate the 
absolute uncertainty estimates due to Wagener et al. 



wovelength, /xm 


Fig. 11.— Aerosol-free Neptune spectra computed for 500 km- 
amagats of H 2 , without He, excluding polarization, and assuming 
Ford and Browne (1973) Raman cross sections (dot-dash for A 5 =l, 
dashed for A5=0) and 2/3 Ford cross sections (solid), compared 
to earlier calculations by Cochran and Trafton (1978) (heavy dot- 
dash) and by Savage et al. (1980) (heavy dashed). Our spectra 
have been smoothed to an effective resolution of 5 nm to facilitate 
comparisons. See text for discussion. 

(1986). The lUE results are 5% to 10% below the theo¬ 
retical clear-sky calculated spectrum, generally showing 
the smaller difference at shorter wavelengths (0.21-0.23 
^m) and the best defined and larger differences at longer 
wavelengths (0.29-0.33 p.m), which is likely due to an 
absorbing high altitude haze. Our Haze I model calcula¬ 
tions match the lUE observations generally well within 
the lUE error bars. If we had neglected polarization, 
the model calculations would drop about 5%, making 
the clear-atmosphere model the best fit near 0.25 ^m. 
In this spectral range, polarization makes the difference 
between needing haze absorption and not needing it. 

The Faint Object Spectrograph (FOS) observed Nep¬ 
tune at a phase angle of 1.22° and a sub-earth lati¬ 
tude of -24.78° planetographic on 19 August 1992 us¬ 
ing apertures of 0.3", covering 21-36°S, and 1.0", cov¬ 
ering 3-35°S. Because these apertures are smaller than 
Neptune’s 2.4" disk and placed near the middle of Nep¬ 
tune’s disk, the observed I/F and the relative ampli¬ 
tudes of Raman spectral features are increased relative to 
what would be observed for a full-disk average (shown in 
Sec. [7]). Courtin (1999) obtained uncalibrated geometric 
albedo spectra by dividing FOS flux spectra by a solar 
flux spectrum measured by the UARS SOLSTICE in¬ 
strument (Rottman et al. 1993; Woods et al. 1993). The 
FOS results we use are the 1-arcsec spectra degraded to 
a resolution of 1 nm to improve the signal/noise ratio 
(shown in Courtin’s Fig. 2). Courtin performed a ra¬ 
diometric calibration by matching his 0.22-0.33 p,m FOS 
spectra to full-disk 0.3 - 1.0 p,m groundbased spectrum of 
Karkoschka (1994), obtained during 23-26 July 1993 at a 
phase angle of 0.4°. The FOS observations thus have the 
same 4% absolute uncertainty as the Karkoschka (1994) 
reference spectrum. The FOS noise level varies from 
0.4% at 0.33 p,m to 1% at 0.26 /rm. Courtin’s Neptune 
spectrum is compared to our clear-sky and Haze I calcu¬ 
lations in Fig. 1131 Our calculations as a function of view 
angle were integrated over approximately the same field 
of view as the FOS observations to account for the in¬ 
creased Raman effect near normal view angles. We then 
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convolved them to a final resolution of 1 nm and scaled 
this average spectrum to match our calculated geomet¬ 
ric albedo spectrum in the wavelength range from 0.31 
fim to 0.33 ^m, to simulate the calibration procedure 
used by Courtin. Note the excellent agreement between 
the shapes and amplitudes of the calculated spectral fea¬ 
tures for the Haze 1 model and those measured by the 
FOS. The most glaring exceptions (a and b in the fig¬ 
ure) appear to be artifacts in the FOS spectrum. These 
anomalies are not present in the overlapping portion of 
the Karkoschka spectrum, which is in good agreement 
with the calculated spectrum. The spikes c and d in the 
ratio spectrum are at points where the 1/F spectrum has 
very sharp gradients, and could be removed by a slight 
wavelength shift. The FOS/Haze 1 ratio spectrum gen¬ 
erally agrees with the lUE/Haze 1 ratio shown in Fig. 
1121 indicating that haze absorption may be depressing 
the reflectivity by about 6 %- 10 % below the theoretical 
value for a clear Neptune atmosphere in the 0.23-0.27 
fim range, with increasing absorption indicated at longer 
wavelengths. 



0.22 0.24 0.26 0.28 0.30 0.32 

wavelength, fj.m 

Fig. 13.— Comparison of our aerosol-free (dot-dash) and Haze I 
(solid) calculated Neptune spectra with FOS observations (-1—|—h) 
by Courtin (1999). All spectra are convolved to a spectral resolu¬ 
tion of 1-nm. 
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Fig. 12.— Comparison of our aerosol-free (solid) and Haze I 
(dashed) calculated Neptune spectra with observations by the lUE 
(dot-dash). The model spectra are smoothed to a resolution of 
5-nm, to be compatible with the lUE observations. Ratio spectra 
are displayed in the lower panel. 


Figure[T4]compares clear-atmosphere and Haze I model 
calculations with the disk-integrated groundbased spec¬ 
tra of Karkoschka (1994). Between 0.35 ^m and 0.45 
fj,m the observations are about 13% below the clear- 
atmosphere calculations, indicating a presence of an ab¬ 
sorbing haze that is even more influential than at shorter 
wavelengths. This behavior suggests that the haze ex¬ 
tends deep enough into the atmosphere that it becomes 
more influential as the molecular scattering optical depth 
above it decreases with increasing wavelength. This 
is qualitatively consistent with microphysical models of 
Moses et al. (1995), which suggest that haze opacity in¬ 
creases down to the 870 mb level. It is also roughly 
the character of our Haze I model, which is shown to 


provide good agreement with Karkoschka’s observations, 
although the computed amplitudes of the Raman spec¬ 
tral features are somewhat larger than the observations 
in several cases. This discrepancy might mean that more 
aerosol opacity is needed in the 1-3 bar level, or that 
there is a slight difference in the effective spectral resolu¬ 
tion of the observations and the calculations. At longer 
wavelengths, the high observed albedo values in the deep 
CH 4 bands at 0.8 /im and 0.9 /rm, indicate that there is 
also particulate scattering at high altitudes, which is not 
included in either of the haze models. Much of that con¬ 
tribution probably comes from discrete bright features 
seen in Voyager, HST, and groundbased images. Effec¬ 
tive pressure estimates for such features range from 23 
mb to 60 mb at 30-40° N, 100-230 mb at 30-50° S, and 
170 mb to 270 mb near 70° S (Sromovsky et al. 2001b; 
Gibbard et al. 2003). The effective fractional coverage of 
bright high altitude clouds required to explain the albedo 
in deep CH 4 bands appears to be about 0.5% to 1% for 
an upper cloud at 150 mb (Sromovsky et al. 2001a). The 
effective fractional coverage of a coexisting 1.3-bar cloud 
is ~1% assuming a unit-albedo Lambertian reflector. 

The observed albedo in the window regions at 0.825 
fim and 0.935 fj.m is also significantly above the clear- 
sky calculation, suggesting a significant aerosol contribu¬ 
tion that could be at mid to deep levels, which is largely 
satisfied by the 1.3-bar and 3.8-bar clouds in the Haze 
I model. The excess I/F calculated in the window re¬ 
gions at 0.59 /im, 0.63 /xm, 0.68 /xm, and 0.75 /xm, can 
be substantially reduced by increasing the CH 4 contin¬ 
uum absorption by amounts comparable to Karkoschka’s 
stated continuum uncertainty. It could also be reduced 
by decreasing the single-scattering albedo of the cloud 
aerosols. The observations seem to show much less of 
the collision-induced absorption in the 0.8-0.83 /xm re¬ 
gion than is evident in the calculations. Another way to 
describe this difference is to say that there is an imbal¬ 
ance in the 0.83 and 0.94 /xm window regions that is evi¬ 
dent in the calculation, but not in the observations. That 
imbalance may be due to an error in the CH 4 absorption 
coefficients. If we add to Karkoschka’s standard coeffi¬ 
cients the continuum absorption uncertainty, we actually 
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obtain very similar albedo values in these two windows. thermal equilibrium. 



Fig. 14. — Comparison of calculated aerosol-free (red) and Haze 
I (green) model geometric albedo spectra of Neptune with ground- 
based observations (black line) by Karkoschka (1994). The ratios 
of model calculations to observed are shown in the lower panel. 


5.6. Spectral Evidence of the Ortho/Para Ratio 

Because the J = 0 Raman cross sections are larger 
than the J = 1 cross sections (Fig. [2]), the ortho/para 
ratio will have an effect on the Raman features observed 
in Neptune’s reflection spectrum. The largest difference 
in ortho/para ratios for equilibrium and normal H 2 oc¬ 
cur near 100 mb (Fig. [T|), and thus associated spectral 
variations are more likely at short wavelengths where 
this region is near the t = 1 level (Fig. [d]). Calcula¬ 
tions (Fig. [T 5 I) show that the largest difference is for the 
peaks at 0.280 /im and 0.2852 /rm, which are about 5% 
larger for equilibrium H 2 . The overlay of FOS obser¬ 
vations suggests better agreement with equilibrium H 2 , 
as concluded by Courtin (1999), although the relative 
size of the smaller peak is very dependent on the effec¬ 
tive resolution of the observations and both peaks have 
amplitudes that depend on view angle. And since the 
FOS observations cover only the central disk, it is not 
strictly valid to compare them to geometric albedo obser¬ 
vations. Uncertainties in the wavelength scale and effec¬ 
tive resolution of the FOS observations arise from com¬ 
parisons with overlapping groundbased observations of 
Karkoschka (1994), which show significantly less spectral 
modulation in the 0.3-0.33 ^m overlap region, as noted 
earlier. The presence of an absorbing haze and high- 
altitude discrete clouds also affects these peaks. Thus 
it seems premature to make a strong conclusion based 
on these observations. However, Courtin’s (1999) result 
that /eg=0.88±0.23 is based on higher resolution obser¬ 
vations than those shown in Fig. 1151 and was derived 
after special processing to minimize wavelength errors. 
Courtin’s results are also compatible with independent 
conclusions by Baines and Smith (1990) and Conrath et 
al. (1991) that the ortho/para ratio for Neptune is near 




Fig. 15. — Geometric albedo spectra for equilibrium H 2 (solid) 
and normal H 2 (dot-dash) compared to the observed FOS spectrum 
(plus signs) (Courtin 1999) scaled by the factor 1.08. The ratio of 
equilibrium to normal albedos is in the lower panel. 


6 . APPROXIMATIONS OF RAMAN SCATTERING 
6 .1. The Karkoschka Correction 

Karkoschka (1994) presented a method for transform¬ 
ing spectra between Raman and non-Raman forms based 
on the assumption that the measured geometric albedo 
spectrum p{v) can be well approximated by a linear com¬ 
bination of terms involving offset versions of the “phys¬ 
ical” spectrum q{v), which is the reflectivity spectrum 
without Raman scattering. The mathematical expres¬ 
sion of this method is 

p{v) = h{v)q{i') + 

T, + ^v)q{v P Ai/)F0(i/ -F 1^v)/Fq{u) (31) 

Av 

where Fq is here the solar photon spectrum, fo^v) is 
the fraction of photons not Raman scattered, /aj/ is the 
fraction undergoing transition Ai/, and where fo{v) + 
"^fAviv) = 1. Equation [31] can be inverted to obtain 
the Raman-free “physical” spectrum from the observed 
spectrum using qi{v) = p(v) as a first guess, and then 
improving the guess using the iterative equation 

qn+iiy) = fa{v)~^[p{v) - 

T, + Azz)g„(i/ -F Av)Fq{v + Ai/)/F0(i/)] ,(32) 

Av 

for which four iterations are usually sufficient to converge 
on q for a given set of fractions. Karkoschka assumed 
a power law A-dependence for /a^, then used minimum 
roughness of the fitted spectrum for 0.31 pm < A < 0.405 
/rm as the criterion for picking the optimum fractions and 
exponent. 

We first investigated how well the proposed model 
could fit Raman spectral calculations and how the fit¬ 
ted constants related to Raman scattering cross sections. 
We considered both a clear-atmosphere model and a 
model with haze and cloud contributions. We first fit 
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Fig. 16. — Application of Karkoschka’s method of removing effects of Raman scattering from the model spectrum (thin solid line) for our 
haze/cloud aerosol structure. The inverted “physical” spectrum (thick solid line) is slightly above the true Raman-free spectrum (dot-dash 
line). The smoothing length used in judging fit quality was 0.71 nm. Lower: the ratio of the fit to the smoothed fit, with the heavier 
portion of the curve indicating the part used to constrain the fit. 


TABLE 5 

Direct Fits of Karkoschka’s Empirical Model to Sample Spectra. 


Spectrum 

A-range 

(pm) 

/o 

/so 

/si 

fa 

n\ 

RMS 

DEV 

Clear Sky 

0.31-.405 

0.100±0.005 

0.065±0.007 

0.050±0.005 

0.085±0.005 

-0.80±0.05 

0.93% 

Clear Sky 

0.23-.405 

0.100±0.005 

0.070±0.005 

0.030±0.005 

0.091±0.005 

-0.50±0.07 

1.47% 

Haze II 

0.31-.405 

0.050±0.005 

0.065±0.005 

0.030±0.005 

0.057±0.005 

-0.55±0.05 

0.44% 

Haze H 

0.23-.405 

0.035±0.003 

0.060±0.004 

0.030±0.005 

0.085±0.005 

-0.35±0.06 

0.90% 

Haze H (inv) 

0.23-.405 

0.060±0.020 

0.060±0.004 

0.030±0.003 

0.065±0.005 

-0.35±0.05 

0.58% 

K94 Coeff. 

0.31-.405 

0.035±0.015 

0.060±0.01 

0.025±0.015 

0.050±0.015 

-1±1 



NOTE: Fractional values are given at 400 nm. Values at other wavelengths are proportional to A">'. Haze II 
(inv) refers to the inversion of fit parameters from the Raman spectrum only, rather than starting with the 
Raman-free spectrum. In this case RMS DEV is relative to the smoothed spectrum, which was smoothed 
with a boxcar of 0.71 nm. The last row is from Karkoschka (1994). 


the fractional parameters and wavelength dependence ex¬ 
ponent by minimizing the RMS deviation between the 
Karkoschka model and the calculated spectrum, first over 
the 0.31-0.405 ^m range used by Karkoschka and also 
over the wider 0.23-0.405 /im range. The results are 
summarized in Table 5. Note that our fit results are 
comparable to those of Karkoschka, especially when ap¬ 
plied to a hazy atmosphere. Although we were able to 
obtain good fits over the limited spectral ranges (RMS 
deviations of 0.44% and 0.9% for the haze/cloud model, 
and 0.93% and 1.47% for clear sky model), the fits did 
not accurately match the observed spectrum at slightly 
longer wavelengths, especially where CHq absorption was 
present. The weak CHq bands in the Karkoschka conver¬ 
sion of the physical spectrum were much stronger than 
the corresponding bands in the correct Raman calcula¬ 
tion. This is understandable as a consequence of the 


nonlinear relationship between 1/F and single-scattering 
albedo. A small absorption added to a conservative at¬ 
mosphere can have a much larger fractional effect than 
it does when added to an absorbing atmosphere. For ex¬ 
ample, adding kabs = 0.003 x kscat to a Rayleigh atmo¬ 
sphere with kabs = 0 changes oj from 1.0 to 0.997, which 
decreases the I/F from 0.791 to 0.699, a 9% drop. But 
if the 1/F is already at 0.559 due to Raman absorption, 
adding the same CH 4 absorption as before would change 
OJ from 0.975 to 0.972, producing only a 2% drop in 1/F. 
Conhrming Karkoschka’s fit results, we find that the vi¬ 
brational coefficient (/q) is about twice the size of the co¬ 
efficient for the J = 1 rotational contribution (/si), even 
though its Raman cross section is somewhat smaller (Ta¬ 
ble 3). Since we use these cross sections to compute the 
spectra, the Ht results do not suggest that there is some¬ 
thing wrong with the cross sections. A large part of this 
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result comes from the fact that the vibrational coefficient 
represents the average of ortho and para contributions, 
while the individual ^(O) and S'(l) contributions are each 
multiplied by fractions that decrease their relative con¬ 
tributions by roughly a factor of two (for fpara ~ 1/2, 
see Eg. [H and Fig.[T]). Additional factors are discussed in 
Sec. 15(11 

We next consider how well Karkoschka’s method is 
able to remove Raman scattering effects from a spec¬ 
trum that includes Raman scattering. We inverted our 
model Raman calculations using the inversion method 
based on Eg. |32l The retrieved parameter values for the 
Haze II test case are given in Table 5. As illustrated 
in Fig. [m we found pretty fair agreement between the 
true Raman-free spectrum and inverted spectrum for the 
haze/cloud case when we used the 0.23-0.405 /rm range 
for the fitting, although our inverted spectrum is offset 
about 0.01 albedo units above the true spectrum, and 
the inverted parameters differ somewhat from the best- 
fit parameters. Using the 0.23-0.405 /rm spectral range 
yielded a smoother inverted spectrum but a larger offset 
albedo of 0.03. Much worse results were obtained with 
the clear sky model, which has larger Raman features 
that are not as well fit by Karkoschka’s empirical model. 
As expected, we found that the features due to weak 
CH 4 bands were much smaller in the transformed spec¬ 
tra than in the true physical spectra. Thus, estimates 
of CH 4 absorption coefficients using spectra corrected in 
this way will be significantly below the true absorption 
level. 

Karkoschka’s method has several physical flaws that 
limit its applicability. Multiple scattering is one effect 
tha t is not accounted for but can be significant (see Sec. 
EU). Another problem with the physics assumed in Eg. 
m is that the contribution at the present wavenumber is 
not necessarily proportional to the geometric albedo at 
the upward shifted wavenumber. In fact, the proportion¬ 
ality constant also depends on the vertical distribution 
of the scattered photons in relation to the vertical dis¬ 
tribution of opacity at the scattered wavelength. The 
Karkoschka forward method also is not generally appli¬ 
cable to spatially resolved observations because the am¬ 
plitude of Raman effects depends on view angle as well as 
local aerosol structure. The useful spectral range for the 
Karkoschka correction is also rather limited, probably to 
wavelengths less than about 0.405 /rm, partly because it 
does not properly modify weak CH 4 absorption bands. 

6 . 2 . The Wallace Approximation 

The approximation suggested by Wallace (1972) is a 
modification of the molecular single scattering albedo of 
the atmosphere by treating rotational transitions as extra 
sources of scattering and the vibrational transition as a 
pure absorption. This can be expressed in terms of cross 
sections as follows: 

^Ray “t” f^rot “t” 

ujw = -f-----, (33) 

where is the Rayleigh scattering cross section, Uabs 
is the total cross section for absorption by CH 4 and CIA 
by H 2 , (Trot is the combined cross section for rotational 
Raman transitions, dvib is the cross section for vibra¬ 
tional Raman transitions, and where /? = 0 for the orig¬ 


inal Wallace approximation. For a clear Neptune atmo¬ 
sphere, we find that the original Wallace approximation 
provides a crude match to the low resolution baseline of 
the spectrum, but is generally 5% to 10% low and of 
course does not produce characteristic Raman spectral 
features. A few detailed calculations by Wallace (1972) 
at 0.2 pm and 0.4 pm also indicated that his approxi¬ 
mation was about 5% low. The Wallace approximation 
error is nearly zero in the near-IR CH 4 windows, but 
is about 4% low in the CH 4 absorption bands, which is 
about what would be expected if Raman scattering were 
ignored (see lower part of Fig. [T^). Figure [HI displays a 
modified form of the Wallace approximation that boosts 
the I/F value by setting /3 = 0.433, which treats the vi¬ 
brational cross section as 56.7% absorption and 43.3% 
scattering. This improves the overall agreement in the 
blue to orange part of the spectrum, but has little ef¬ 
fect elsewhere. The Wallace approximation approaches 
the exact result when high clouds obscure most of the 
molecular scattering. 
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Fig. 17.— Upper: Geometric albedo spectrum computed for a 
clear Neptune atmosphere with no accounting for Raman scattering 
(blue), with Raman scattering included (black), and with Raman 
scattering approximately accounted for using a modification of the 
Wallace (1972) approximation (violet), a modified and generalized 
version of the Pollack (1986) approximation (blue) and the tuned-oj 
approximation (red). Lower: Ratios of the approximate geometric 
albedo spectra to the accurate Raman calculation. 


6.3. The Pollack Approximation 

Pollack et al. (1986) approximated Raman scattering 
contributions from shorter wavelengths by scaling the 
Raman scattering cross section at the reflected wave¬ 
length by the ratio of solar irradiance at the source wave¬ 
length to that at the reflected wavelength. Generalizing 
this approach to include more than one transition, the 
molecular single-scattering albedo in this approximation 
can be written as: 

U}p[v) = min ---i— -=-, 1 

CTRayW) + cra,hsW)+2^ecrKeW) > 

(34) 

where ri = FqIv -|- /S.V()jF q{v), gi is a function of rg 
that can be empirically adjusted to compensate for some 
of the complexities that are glossed over by this con¬ 
cept, and the upper bound of 1 is used to avoid energy 
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conservation problems. Setting gvib(''’vib) = ?'vib, and 
omitting the rotational transitions completely, would re¬ 
duce this to Pollack’s original suggestion, while setting 
( 7 rot = 1 and ^vib = P would reduce this to the Wallace 
approximation. A problem with Pollack’s original for¬ 
mulation was that it resulted in wp > 1 at some wave¬ 
lengths, as noted by Courtin (1999), which leads to in¬ 
stability in the solution of the radiation transfer equa¬ 
tion. Another problem is that it considered only the Qi 
vibrational transition, which alone would be unable to 
generate the ghost features arising from the rotational 
transitions. The form given in Eg. 1341 has the capability 
to fit all the main ghost features and has the potential, 
with tuning of the form of ri, to yield the correct baseline 
value as well. It also has an advantage over Karkoschka’s 
correction, in that it is introduced at a more fundamental 
physical level that is naturally modified by the presence 
of aerosols. However, wherever the geometric albedo ex¬ 
ceeds the value for a conservative Rayleigh atmosphere 
of 0.7908 (Sromovsky 2004), even ujp = 1 will not be 
able to reproduce the albedo value. This problem can 
be reduced by working at lower spectral resolutions for 
which these extreme peak amplitudes become reduced 
below 0.7908. 

Preliminary calculations with this approximation used 
the three transitions given in Table 4, with gg = ciri and 
ce the same for all transitions. We found that making 
the constant Cf = 0.875 for all i produced a reasonably 
good match to the baseline level of geometric albedo, but 
generated spectral modulations that were too large. To 
limit these modulations, we used an exponential function 
to limit the effect of the spectral irradiance ratio, i.e. we 
used gi = w(l — exp(—r^/rc)), which is ~ rg for small ri 
and approaches w for large values of r^. We found that 
w = 2 provided a reasonable fit to the Raman spectral 
features and to the baseline level, as shown in Fig. flTl 
While this is certainly an improvement over the Wallace 
approximation, and does reproduce many of the spec¬ 
tral features associated with Raman scattering, it does 
not do so consistently and accurately enough to use for 
interpreting those features. Like the Wallace approxima¬ 
tion, it also fails to improve geometric albedo accuracy 
for A > 0.47 /rm. Both Wallace and Pollack modifications 
to the single-scattering albedo seem to be far less effec¬ 
tive in regions where CH 4 absorption is important. Since 
the actual source function depends on the light level in 
the atmosphere, which may have a spectral shape that 
is very different from the incident solar spectrum, it is 
not appropriate to use the solar spectral ratio as a scale 
factor in spectral regions with strong CH 4 absorption. 
The following method goes part of the way in solving 
this problem. 

6.4. Approximation of Raman Scattering by Spectral 
Tuning of uj 

The following describes a new approximation of Ra¬ 
man scattering that uses a different kind of modifica¬ 
tion of the molecular single scattering albedo. Instead 
of using a constant factor or the solar spectral ratio to 
modify the molecular single-scattering albedo, this ap¬ 
proximation uses a A-dependent multiplier that is tuned 
to minimize the difference between simulated and exact 
Raman calculations of geometric albedo for an appropri¬ 
ate model atmosphere. The fundamental relation used is 


from Sromovsky (2004): 


p{uj) = 0.7908a;-°-2®®[l-exp(-|logio(l-a;)|i-2®Vl-504)], 

( 35 ) 

where p is the geometric albedo including polarization 
contributions and uj is the single scattering albedo of a 
semi-infinite Rayleigh atmosphere. This equation is in¬ 
verted numerically to provide a function uj{p) so that a 
given geometric albedo can be related to an equivalent 
single-scattering albedo of a semi-infinite Rayleigh at¬ 
mosphere. Starting with a Raman-free calculation pro¬ 
ducing a spectrum pffpiv), and a true Raman calcu¬ 
lation producing a spectrum pp^u), we use the inverse 
transform to compute equivalent single-scattering albe¬ 
dos ujp^NR and ujp^R. Noting that l/w = 1-1- kabs/kscatt, 
we then attribute the difference between the equivalent 
inverse single scattering albedos as due to a difference in 
the effective absorption due to Raman scattering. Thus 
we are led to the approximation for modifying the sin¬ 
gle scattering albedo of the molecular scatterings that 
simulate true Raman scattering, namely, 


UJ 


sim,n-\-l — 


l/(^szm,n 


(36) 


where ujsim,n is the single-scattering albedo simulation of 
Raman scattering on the prior iteration, ujp^n is the effec¬ 
tive semi-infinite single scattering albedo obtained from 
the spectrum generated using uJsim,n for the molecular 
scattering. Because Neptune’s atmosphere has a ver¬ 
tically inhomogeneous CH 4 distribution, and thus will 
not satisfy Eq. |35l it is necessary to iterate a few 
times to achieve convergence. On the first iteration, 
we have ujsim,n=i = 1-0 and a;p,„=i = ujp^NR- The 
molecular single-scattering albedo spectrum that results 

= ani7R)ZZsH ’ reduces to ujt = uj,im 

when (Tabs = 0 . 

The nearly perfect agreement between the geometric 
albedo computed with Raman scattering and the approx¬ 
imation using a tuned single-scattering albedo is shown 
in Fig. 1171 Unlike the previously discussed approxima¬ 
tions, it does provide useful accuracy for A > 0.45 /im. 
Like the other approximations, it cannot reproduce the 
peaks that exceed the I/F for a unit geometric albedo. 
While the modified single-scattering albedo can exceed 
unity, or even go negative, that does not lead to prob¬ 
lems as long as the other contributors to the scattering 
process lead to a combined single-scattering albedo < 1 . 

The important question is whether this approximation 
or any of the others that are similarly constructed, are 
useful with clouds present, and for center-to-limb scans. 
Because all of these approximations that modify single¬ 
scattering albedos really lack the essential physics of the 
Raman scattering process, it is doubtful whether they 
can have much utility in situations for which they have 
not been tuned. A partial test of this utility is the degree 
to which they can simulate the angular variation of Ra¬ 
man scattering, which is described in the next section. 
Section | 8 ] compares their performances when haze and 
cloud aerosols are present. 

7. ANGULAR DEPENDENCE OF RAMAN SCATTERING 
EFFECTS 

Raman scattering effects on Neptune’s I/F spectrum at 
zero phase angle are strongest near the center of the disk 
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(zero zenith angle) and weakest near the limb fFig. [T51) . 
The errors obtained by ignoring Raman scattering (not 
shown) are quite substantial at moderate view angles, 
especially at short wavelengths and at high spectral res¬ 
olution, ranging from about 20 to 60% in the 0.25-0.3 p,m 
range at a zenith angle of 8.1°. At view angles of ^60°, 
the fractional errors are comparable to those seen for the 
geometric albedo. As noted for the geometric albedo, 
errors ^4% are seen in the deep CH 4 bands, with very 
little error seen in the near-IR window regions. In the 
longer wavelength bands, the fractional error does not 
improve much near the limb, unlike the error at shorter 
wavelengths, which decreases substantially. 

The modified Pollack approximation helps consider¬ 
ably in reducing errors at short wavelengths (Fig. [TSll . 
The extreme error near the central disk drops from 62% 
to 14% and the RMS error from 19% to 4%. At 63° 
the RMS error drops from 8 . 8 % to 2.2%. But at long 
wavelengths, the Pollack approximation actually makes 
things slightly worse in the middle of strong absorption 
bands, probably because the vertical location of Raman 
photons becomes more critical to the resulting I/F level. 
Overall, the modified Pollack approximation has value 
in modeling a variety of observations. It is far better 
than ignoring Raman scattering and is roughly twice as 
accurate as the modified Wallace approximation. How¬ 
ever, using the modified Pollack approximation can lead 
to substantial errors in limb-darkening profiles, depend¬ 
ing on the wavelength that is considered. Sample profiles 
are shown in Fig. [T21 

Under most conditions, the tuned-w approximation is 
even more accurate than the modified Pollack approx¬ 
imation, as evident from comparison of spectra at the 
three sample angles (Fig. fT^ and from a comparison of 
center-to-limb scans at sample wavelengths (Fig. IT^ . 
The extreme error near the central disk drops to 10.3% 
and the mean and RMS deviations to only -1.1% and 
1.2%, with comparable results at 63° view angles. The 
worst errors are seen at large view angles. At 85.7° , the 
extreme error increases to 24% and the mean and RMS 
errors increase to 2.6% and 3.3%. But most of these 
errors occur in the CH 4 bands beyond about 0.53 ^m. 
If restricted to 70° view angles and A < 0.85 /rm, this 
approximation is remarkably effective, with errors rarely 
exceeding a few percent. 

8 . EFFECTS OF AEROSOLS ON APPROXIMATION 
PERFORMANCE 

Calculations for the Haze H atmospheric model con¬ 
taining the cloud and haze aerosols described in Sec. 14.21 
are presented in Fig. [20l where the calculation that ig¬ 
nores Raman scattering is compared to the Raman calcu¬ 
lation and to the modified Pollack and tuned-w approxi¬ 
mations. At UV-visible wavelengths we see that the size 
of the Raman features is substantially reduced by the 
addition of a high altitude absorbing haze, which also 
depresses the baseline reflectivity as needed to approx¬ 
imately match observations. The net effect of Raman 
scattering on the baseline reflectivity is also reduced sub¬ 
stantially compared to its effect in a clear atmosphere. 
This is another example of the nonlinear effects of adding 
absorption (see discussion in Sec. OD . The effects of the 
cloud layers are most easily seen at longer wavelengths 
where the I/F in window regions and regions of inter- 





Fig. 18.— Comparison of Neptune’s I/F spectra computed with 
Raman scattering (black) and with a modified Pollack approxima¬ 
tion (blue) and tuned-oj approximation (red), for three observing 
angles at zero phase. The ratios of the two spectral versions are 
shown at the upper part of each panel; in the top panel the ratio 
is offset by 0.4 to prevent overlap. The spectra are convolved to a 
resolution of 36 cm“^. 

mediate methane absorption are increased substantially 
compared to that computed for a clear atmosphere. 

The tuned-w approximation is seen to perform well for 
A < 0.6 /im, where haze effects dominate, but is biased 
about 5% too high in the strongly absorbing regions at 
longer wavelengths, although remaining quite accurate 
in window regions. This is not too surprising, because 
the tuning in these regions is involves a delicate balance 
between large molecular absorption and large corrections 
of that absorption to simulate the Raman-scattered con¬ 
tributions. As vertical structure is changed from that 
used for tuning, the tuning can be easily upset. 

The modified Pollack approximation performs some¬ 
what worse at UV-visible wavelengths than it did for 
the clear atmosphere case (Fig. [T7]). The overall nega¬ 
tive bias seen here could probably be largely eliminated 
by adjusting the parameters of the gt, factors in Eq. [M) 
That would do much to reduce the positive errors in the 
window regions beyond 0.45 /im or the negative errors in 
the absorbing regions. These are inherent in the physi¬ 
cal model on which the Pollack approximation is based. 
Overall, this approximation is less upset by the presence 
of the cloud layers than is the tuned-w approximation. 

A comparison of these two approximations at three dif¬ 
ferent view angles is presented in Fig. We see that 
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Fig. 19.— Comparison of Neptune’s center-to-limb profiles at 5 
selected wavelengths computed with Raman scattering (solid) with 
a modified Pollack approximation (blue) and tuned-oj approxima¬ 
tion (red). The ratios of the two spectral versions are shown in the 
lower panel. The calculations were convolved to a resolution of 36 
cm“^ before sampling the selected wavelengths. 


the tuned-w approximation is much better near the cen¬ 
tral disk, and effective at all wavelengths. It is also ef¬ 
fective at all angles for wavelengths less than 0.47 ^m. 
But at longer wavelengths, errors in the regions of strong 
methane absorption grow with increasing view angle, 
making the modified Pollack approximation better for 
view angles beyond about 50° . 

The Karkoschka approximation was not included in 
this particular comparison because is not useful for spa¬ 
tially resolved observations, especially those with vari¬ 
able cloud structure, unless the model coefficients are 
changed for each view angle and structure. Even then, 
the serious problems at A > 0.4 ^m (see Fig. [TH]) also 
reduce its utility, even for disk-average analyses. 

9. SUMMARY AND CONCLUSIONS 

Major results of this investigation can be summarized 
as follows: 

(1) Radiation Transfer.: Procedures for modifying the 
Evans and Stevens (1991) vector radiation transfer 
code to handle Raman scattering are described and 
justified. Equations are provided for the Raman 
differential generators. The use of a Raman source 
matrix facilitates the computation of the first spa¬ 
tially resolved Raman scattering spectra for Nep¬ 
tune. Validation of the code is obtained by demon¬ 
stration of photon conservation, comparison with 
prior low-resolution calculations, and comparison 
with HST and groundbased observations. 
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Fig. 20.— Geometric albedo of Neptune for a sample aerosol 
structure containing an absorbing high-a ltitu de haze and two lay¬ 
ers of clouds (Haze II model, given in Sec. l4.2( t. computed with Ra¬ 
man scattering ignored (green), using the Pollack approximation 
(blue), the tuned-o; approximation (red), and the exact method of 
calculation (black). The ratios to the true calculation are displayed 
in the lower panel. 
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Fig. 21.— Comparisons of Neptune’s spectra at specific view an¬ 
gles for a sample aerosol structure containing an absorbing high- 
altitude haze and two layers of clouds (Haze H model), using the 
modified Pollack approximation (blue), the tuned-o; approxima¬ 
tion (red), and the exact method of calculation (black). The ratio 
spectral are also shown in each panel. 


(2) Basic Effects on Geometric Albedo.: Raman scatter- 
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ing reduces the baseline geometric albedo of a 
semi-infinite conservative Rayleigh scattering at¬ 
mosphere by about 25% from the value of 0.791 
it would otherwise have. In the presence of a high- 
altitude absorbing haze, the effect of Raman scat¬ 
tering is reduced by about a factor of two. The 
well known sharp positive spectral peaks from fill- 
in of absorption lines in the solar spectrum can 
greatly exceed the Rayleigh conservative limit. For 
photons introduced at 228 nm, the Q transition 
source function for the first scattering has a broad 
peak near 300 mb; after further scattering photons 
move deeper into the atmosphere and leak out the 
top, moving the peak of the source function further 
downward after each scattering. 

(3) Near-IR Effects.: For a clear atmosphere about 4% 

of the light in the deep methane bands in the near- 
IR is due to Raman scattering, while the near IR 
window regions are almost unaffected by Raman 
scattering. The conclusion of Cochran and Trafton 
(1978) that the light in the deep methane bands 
is almost entirely Raman scattered is a result of 
their assumption of what are now known to be ex¬ 
cessively high methane mixing ratios in the strato¬ 
sphere. 

(4) Multiple Scattering.: The calculated history of 

monochromatic incident photons shows that 
multiple Raman scattering is quite significant in 
a conservative Rayleigh scattering atmosphere, 
where the Raman absorption process can compete 
with the loss to space. In one example calculation, 
only 40% of the Raman scattered photons exiting 
the atmosphere experienced only one Raman 
scattering. 

(5) Angular Variation.: At low phase angles spatially re¬ 

solved Raman spectra show that Raman spectral 
features are enhanced at near the center of the 
planetary disk and suppressed at the limb. We also 
found that approximations that work reasonably 
well in matching accurate disk-integrated calcula¬ 
tions do not do as well in matching observations in 
specihc directions, especially at large view angles. 

( 6 ) Comparison with Observations.: Aerosol-free models 

of Neptune’s spectrum correlate well with observed 
spectral features, but confirm a need for haze ab¬ 
sorption to reduce the baseline geometric albedo 
and to reduce the amplitude of spectral modu¬ 
lations induced by Raman scattering. lUE full- 
disk observations, adjusted to account for current 
knowledge of of Neptune’s size, have a baseline off¬ 
set of ~ 6 % in the 0.24-0.26 /rm range, increasing 
to about 10% in the 0.30-0.32 ^m range, both of 
which are fit well by our Haze I model, containing 
0.2 optical depths (at 0.5 ^m) of a UV absorbing 
haze of 0 . 2 -/im radius particles uniformly mixed be¬ 
tween 0.1 and 0.8 bars. FOS observations between 
0.22 and 0.29 /xm, have a similar offset and high- 
resolution spectral structure that is also well fit by 
our Haze I model. Groundbased disk-integrated 
observations between 0.35 and 0.45 ^m have a base¬ 
line offset that is 13% below the clear-atmosphere 
calculation which is well matched by our Haze I cal¬ 
culation, although several of the spectral features 


are seen to have greater amplitude in the calcu¬ 
lations than in the observations. We consider the 
Haze I model to be a sample calculation, not a 
tightly constrained fit. 

(7) Karkoschka correction.: The Karkoschka (1994) 

method of applying Raman corrections to calcu¬ 
lated spectra and removing Raman effects from 
observed spectra is shown to be relatively accurate 
under restricted conditions, but has practical lim¬ 
itations for problems involving spatially resolved 
observations or for A > .4 /xm. His Raman re¬ 
moval algorithm undercorrects the depths of weak 
CH 4 absorption bands, and thus corresponding 
absorption coefficients derived from such spectra 
(Karkoschka 1994; 1998) will need some revision. 
The relatively large amplitude of the Q-branch 
contributions found by Karkoschka is shown to be 
consistent with current estimates of Raman cross 
sections, not an indication that they are in error. 

( 8 ) Wallace Approximation.: We agree with Wallace 

(1972) that his suggested approximation yields ge¬ 
ometric albedos that are typically ~5% lower than 
an accurate calculation. A generalization of this 
approximation that counts part of the vibrational 
cross section as contributing to scattering can re¬ 
move this bias. This approximation does not help 
reduce errors in the strong CH 4 bands and is 
mainly useful at low resolution and for A < 0.45 
/xm. 

(9) Pollack Approximation.: The initial form of the ap¬ 

proximation suggested by Pollack et al. (1986) gen¬ 
erates serious errors. But generalizing it to include 
rotational and vibrational transitions, enforcing en¬ 
ergy conservation, and limiting the solar spectral 
ratio enhancement, can produce a relatively good 
approximation, except for A > 0.45 /xm, where ig¬ 
noring Raman scattering is just as accurate. 

(10) Tuned-tu Approximation.: By tuning the absorption 
attributed to Raman scattering it is possible to 
make a non-Raman calculation match a true Ra¬ 
man calculation with virtually no error. If that 
matching is done with respect to geometric albedo, 
it will not match nearly as well for spatially re¬ 
solved calculations. It is nevertheless the best ap¬ 
proximation to use under cloud-free conditions for 
zenith angles less than 70° and wavelengths less 
than 0.85 /xm. It is the only one that provides im¬ 
provements for A > 0.45 /xm, although is not very 
effective for A > 0.47 /xm when cloud layers of in¬ 
termediate opacity are present. 

Most of these conclusions apply with equal force to Ra¬ 
man scattering in the atmosphere of Uranus because of 
its similarity to Neptune in vertical temperature struc¬ 
ture, gas composition, and relatively low aerosol loading. 
Conclusion ( 6 ) would need the most revision to accom¬ 
modate differences in detailed observational results. 
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